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SUMMARY 


This  document  explains  the  need  for  a  new  computer  data  format,  not 
previously  available,  for  the  storage  of  cartographic  information,  v 

cjA  suitable  format  is  described,  together  with  a  set  of  computer  programs 
which  have  been  written  to  enable  data  files  in  this  format  to  be  created,  manipu¬ 
lated  and  displayed  in  various  ways.  A  ^ 
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INTRODUCTION 


Computer  data  can  be  stored  in  a  vatejj^ty  of  ways  and  each  method  is  likely 
to  have  its  own  advantages  and  limitations/^ This  document  is  concerned  with  the 
storage  of  data  representing  some  form  of  'picture'  of  the  Earth's  surface,  as 
viewed  almost  vertically  downwards  from  one  of  the  'Landsat'  series  of  orbiting 
satellites,  using  its  multispectral  scanner, (MSS) . 


The  Landsat  MSS  is  described  in  detail  in  Ref  1.  A  simplified  description 
will  suffice  to  convey  the  information  needed  here.  As  the  satellite  moves  above 
the  Earth,  its  surface  is  scanned  by  means  of  an  oscillating  mirror,  which 
deflects  light,  via  an  optical  system  onto  a  sensor  system.  The  mirror  oscilla¬ 
tion  is  perpendicular  to  the  path  of  the  satellite,  so  that  a  raster-scan  is 
made,  and  during  each  west  to  east  scan  the  output  of  the  sensor  is  digitized  at 
regular  intervals  of  time.  Simplifying  somewhat,  the  result  is  that  a  swath  of 
land  is  imaged,  to  form  a  raster-scan  of  radiometric  values  with  a  constant 
number,  L  of  values  per  scan.  If  the  process  is  stopped  after  M  scans  have 
taken  place,  the  result  is  a  series  of  V  data  values,  where  V  ■  L  x  M  . 


I . I  The  'image'  data  format 

It  is  convenient  to  store  these  V  values  in  the  sequence  in  which  they 
were  obtained,  to  form  a  data  file.  Provided  that  the  number  of  values  per  scan 
is  known  it  is  then  possible  to  reconstruct  the  picture  by  a  suitable  raster-scan 
display,  with  each  picture  element  (pixel)  having  a  brightness  corresponding  to 
its  measured  value.  This  illustrates  a  simple  and  useful  form  of  data  storage  - 
a  sequence  of  values,  which  is  known  to  represent  a  raster-scan.  It  should  be 
noted,  however,  that  the  number  L  of  values  per  scan  must  be  known  (otherwise 
the  picture  cannot  be  reconstructed)  and  preferably  M  also,  so  that  the  size 
of  the  file  (and  picture)  is  known  in  advance.  It  is  therefore  usual  to  have  some 
form  of  information  -  a  'header  record'  -  before  the  main  body  of  the  data,  so 
that  the  latter  can  be  properly  handled. 


A  format  as  described  is  now  in  use  in  Space  Department  of  RAE.  In  detail 
each  data  file  consists  of  a  continuous  sequence  of  8-bit  values.  The  first 
88  such  values  constitute  the  header  information,  and  include  a  value  L  for  the 
number  of  pixels  per  scan  line.  (In  practice  the  number  M  of  scan  lines  is 
also  included.)  The  following  data  consist  of  M  sets  of  L  values,  the  'bright¬ 
ness'  of  the  pixels  of  the  scene.  Such  data  files  are  conventionally  referred  to 
as  image  files  or  images,  and  the  programs  which  handle  them  are  conventionally 
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prefixed  by  the  code  IM  .  Thus,  for  example,  program  IM.LOOK  is  used  to  'look 
at'  an  image  file  in  a  certain  way. 

The  data  storage  method  described  above  has  several  advantages;  in  parti¬ 
cular  it  is  simple  to  use  and  economical  of  storage  space.  However,  it  also  has 
some  limitations,  such  as: 

(a)  Each  pixel  is  located,  by  implication  at  a  specific  place,  an  integer 
number  of  pixels  along  and  down  the  picture.  Fractional  locations  cannot  be 
referred  to. 

(b)  It  is  not  easily  possible  to  establish  any  relationships  between 
pixels,  beyond  their  obvious  sequential  nature.  For  example,  pixels  in  adjacent 
scan  lines  are  not  obviously  related.  Any  feature  on  the  original  scene  which 
extends  over  more  than  one  scan  line  is  not  recognisable  as  such  from  the  data, 
(though  it  may  be  readily  observed  by  human  beings  from  the  resulting  picture 
display) . 

1 .2  The  'map*  data  format 

Considerable  work  is  being  done  worldwide  on  the  interpretation  of  Landsat 

data.  Sometimes  the  resulting  information  can  be  conveniently  stored  in  an 
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image  file.  Some  recent  interpretation  work  done  at  RAE  has  been  concerned 
with  the  extraction  of  cartographic  features,  eg  coastlines,  and  these  are  not 
suitable  for  storage  in  image  files,  in  particular  for  the  two  reasons  mentioned 
in  (a)  and  (b)  above.  Consequently,  it  was  necessary  to  establish  a  more  suitable 
data  storage  format. 

Briefly,  the  new  data  format,  referred  to  as  a  'map'  format,  consists  of  a 
header  record,  followed  by  a  number  of  'features',  each  feature  consisting 
principally  of  either  one  coordinate  value  (to  represent  a  single  location  - 
a  'point'),  or  a  string  of  coordinates  (to  represent  a  continuing  'line').  Such 
a  map  data  format  is  capable  of  storing  much  of  the  diagrammatic  information 
extracted  from  a  Landsat  image. 

This  document  describes  the  map  data  format  in  detail,  together  with  a 
description  of  a  suite  of  computer  programs  and  subroutines  which  has  been 
written  to  handle  such  data  files.  Programs  which  handle  map  files  are  conven¬ 
tionally  prefixed  by  the  code  MAP.,  eg  the  program  MAP. LIST  is  capable  of 
listing  the  coordinates  of  a  map  file. 

The  expression  'line'  or  'linear  feature'  of  a  scene  is  here  used  to  refer 
to  a  sequential  set  of  coordinates  (x,  y  pairs).  If  the  comnencing  and  final 
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coordinates  of  a  line  have  the  same  coordinate  values,  in  both  x  and  y  ,  then 
the  line  is  closed,  ie  it  forms  a  loop.  Because  the  connections  between  each 
coordinate  pair  can  only  be  considered  as  straight  lines,  a  loop  is  in  fact  a 
polygon.  If  the  individual  coordinates  are  sufficiently  close  together,  the 
linear  feature,  when  displayed  on  some  suitable  device,  such  as  an  X-Y  plotter 
or  a  graphics  visual  display  unit,  appears  to  the  eye  as  an  acceptable  representa¬ 
tion  of  a  terrestrial  feature,  that  is  to  say,  the  display  resembles  a  map  etc. 

2  THE  MAP  DATA  FORMAT 

A  map  data  file  consists  of  a  header  record  followed  by  the  body  of  the 
data,  in  sequence.  For  convenience,  and  with  the  objective  of  maximising  the 
reading  and  writing  speed,  the  file  is  divided  into  'blocks',  each  of  1024  words, 
since  that  is  the  number  of  data  words  per  storage  block  of  the  disc  unit  of  the 
host  computer  (Prime)  on  which  the  file  is  to  be  stored.  Should  files  be  written 
to  other  machines,  the  block  size  should  be  changed  suitably,  and  provision  has 
been  made  for  such  a  change  to  be  done  with  the  minimum  of  work. 

The  header  record  is  therefore  1024  words  in  length,  and  the  body  of  data 
is  divided  into  similar  sized  blocks,  although  individual  data  features  may  be 
larger  or  smaller  than  1024  words.  The  Prime  computer  uses  16-bit  words,  and 
the  map  format  has  been  designed  accordingly. 

2. 1  The  header  record 

The  header  record  has  1024  words  available  to  it,  and  this  is  considerably 
in  excess  of  the  amount  normally  needed.  This  space  is  allocated  as  shown  in 
Table  1,  the  details  of  which  are  as  follows. 

2.1.1  'Permanent*  header  description 

Programs  which  create  a  map  file  normally  call  the  subroutine  PERMHD  which 
requests  from  the  user  the  permanent  header  details,  these  being  the  general 
project  name  IPROJ  ,  up  to  12  characters,  plus  the  description  of  the  project, 
IDESCR  ,  up  to  32  characters.  These  details  are  not  changed  by  subsequent  pro¬ 
cessing  of  the  file,  but  are  permanently  available  for  reference. 

2.1.2  Map  projection  details 

Since  the  map  format  is  capable  of  holding  (and  does  on  occasion  hold)  data 
which  is  the  computer  description  of  a  map,  it  is  useful  to  have  the  facility  to 
record  some  information  relating  to  the  cartographic  projection  used. 
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Some  or  all  of  the  following  information  may  be  recorded: 

(a)  IPRJC1:  a  code  number  which  indicates  the  type  of  projection  being 
used.  Codes  so  far  allocated  are: 

IPRJC1  Projection 

1  British  National  Grid 

2  Universal  Transverse  Mercator 

99  Image  file  units. 

The  code  99  is  not  strictly  a  map  projection,  since  it  is  not  cartographically 
defined.  It  is  related  to  the  movement  of  the  Landsat  satellite.  Methods  exist 

for  transforming  'maps'  in  this  pseudo-projection  into  accepted  cartographic 

projections  however,  so  it  is  useful  to  label  a  file  with  the  projection  code  99. 

(b)  IPRJC2:  a  projection  subcode.  This  number  may  be  used  to  qualify 
the  main  code.  The  only  meaning  so  far  allocated  is  for  the  case  when 

IPRJCI  =  2  ,  -ie  UTM.  in  such  cases  it  is  necessary  to  specify  the  zone,  (1  to 

60)  and  this  value  is  held  in  IPRJC2  . 

(c)  SCALE:  this  real  number  is  the  scale  of  the  map.  For  example,  for  a 
one-inch  to  the  mile  map,  SCALE  would  have  the  value  63360.0, 

(d)  XSCALE,  YSCALE:  not  at  present  in  use.  If  it  should  be  necessary  to 
adjust  the  size  of  a  map  to  fit  its  known  corners,  the  X  and  Y  adjustment 
factors  would  be  stored  in  these  locations.  Normally  their  values  would  be  close 
to  unity. 

(e)  SCLLNG,  SCLLAT:  these  numbers  are  available  for  holding  longitude  and 
latitude  values,  usually  stored  as  degrees.  For  example,  in  a  UTM  map,  SCLLNG 
could  contain  the  value  of  the  central  meridian  of  the  zone;  for  a  Mercator  map, 
SCLLAT  could  contain  the  value  of  the  latitude  at  which  the  SCALE  is  accurate. 

(f)  RADEA,  RADEB:  the  Earth's  semi-major  and  semi -minor  axes  in  metres. 
These  values  are  not  at  present  in  use,  but  are  needed  if  map  projections  are  to 
be  considered  in  detail,  or  if  maps  are  to  be  transformed  from  one  projection  to 
another. 

2.1.3  Map  file  statistics 

This  set  of  values  holds  statistics  about  the  map  file,  which  are  useful 
for  various  purposes.  , 

(a)  NLIN:  the  number  of  line  items  in  the  file. 
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(b)  NPT:  the  number  of  point  items  in  the  file.  Since  a  map  file  should 
always  contain  four  corner  points  or  'fiducials',  NPT  should  never  be  less  than  4. 

(c)  NSTEPS,  NSTEPZ:  the  total  number  of  line  ’steps'  in  a  file  is  held 
in  these  two  values.  Any  one  line  may  have  N  coordinate  values,  including 
its  start  and  finish,  and  thus  has  N-l  ’steps'.  The  number  of  steps  for  each 
line  is  totalled,  and  the  value  stored.  Since  this  value  may  be  too  large  to 
hold  in  the  16-bit  integer  available,  it  is  stored  as  two  integers,  modulus 
10000.  That  is  to  say,  if  NTOT  is  the  total  number  of  steps,  then 

NTOT  =  (NSTEPZ  *  10000)  +  NSTEPS  . 

It  would  have  been  possible  to  use  a  REAL  word  or  an  INTEGER*4  to 
store  this  information.  The  chosen  method  allows  larger  numbers  to  be  handled, 
and  more  conveniently,  than  a  REAL  .  INTEGER*4  is  not  available  as  a  standard 
FORTRAN  facility  (though  it  is  becoming  commonly  available). 

(d)  MNFC:  each  item  of  the  map  file  can  carry  up  to  17  'feature  codes', 
which  give  information  about  the  type  of  feature  the  item  represents.  MNFC 
records  the  largest  number  of  such  codes  carried  by  any  item  in  the  file. 

(e)  NDIFFC:  a  file  may  contain  many  items,  and  each  item  may  carry 
several  codes.  The  total  number  of  different  feature  codes  in  a  file  is 
recorded  by  NDIFFC  .  Usually  NDIFFC  is  larger  than  MNFC  . 

2.1.4  Corner  points 

A  map  is  a  representation,  on  say  a  sheet  of  paper,  of  a  portion  of  the 
surface  of  the  Earth.  Two  sets  of  corner  points  exist  therefore:-  the  corners 
of  the  map  and  the  corresponding  locations  on  the  Earth's  surface.  The  eight 
real  words  of  array  CNR(2,4)  hold  the  x  and  y  coordinates  of  the  four 
corners  of  the  map,  reading  from  the  south  west  corner  in  a  clockwise  direction. 

CNR(1,1)  x  coordinate  of  SW  corner 

CNR(2, 1 )  y  coordinate  of  SW  corner 

CNR(1,2)  x  coordinate  of  NW  corner 

CNR(2,4)  y  coordinate  of  SE  corner. 

As  a  convention,  the  SW  corner  is  taken  as  the  origin  of  the  map,  so  its 
coordinates  are  0,  0  .  Also,  the  map  is  defined  as  a  rectangle.  With  these 
restraints,  considerable  redundancy  exists,  and  in  fact  all  of  the  required 
information  may  be  derived  from  the  NE  corner. 


The  CNR  values  are  stored  as  millimetres.  Arbitrarily,  a  maximum  map 
size  of  2000  nm  square  may  be  defined. 

Corresponding  to  the  four  corners  on  the  map  are  four  geographical  locations 
on  the  Earth's  surface.  These  locations  are  stored  in  the  array  GCNR(2,4),  and 
they  correspond  with  the  CNR  values,  ie  GCNR(a,b)  corresponds  with  CNR(a,b)  . 

If  the  map  is  in  British  National  Grid,  the  GCNR  values  are  in  kilometres 
on  that  grid.  For  other  map  projections  GCNR  values  are  held  as  the  relevant 
longitude  and  latitude,  expressed  in  radians.  For  IPRJC1  =99,  ie  the  pseudo¬ 
projection  based  on  image  file  units,  GCNRs  have  no  defined  meaning  at  present. 

2.1.5  Program  logging  information 

The  header  record  has  a  logging  facility  to  allow  a  record  to  be  made  of 
all  processing  done  to  a  map  file.  Some  programs  are  available  to  modify  exist¬ 
ing  files,  and  it  is  useful  to  keep  a  log  of  what  has  been  done.  All  such 
programs,  together  with  the  program  which  initially  created  the  file,  use  a 
subroutine  UPDATE  ,  which  keeps  a  log  in  the  header  record. 

Subroutine  UPDATE  gathers  information,  some  from  the  user,  but  mostly 
from  the  computer  itself  and  from  the  program  in  use,  and  places  this  information 
into  a  40  word  portion  of  an  array  (only  29  of  these  words  are  at  present  in  use). 
There  is  space  in  the  header  record  for  15  such  sets  of  40  words.  After  a  40  word 
set  has  been  assembled,  UPDATE  moves  all  the  preceding  sets  one  place  'down', 
and  then  places  the  latest  set  at  the  top  of  the  available  space.  Thus,  if 
15  sets  of  log  data  had  been  in  use  prior  to  this  operation,  the  oldest  one  would 
have  been  lost,  as  the  others  moved  down.  In  practice,  it  is  not  expected  that 
there  will  be  as  many  as  15  UPDATES  on  most  data  files,  with  the  suite  of  pro¬ 
grams  at  present  available. 

Within  each  40  word  set  of  data,  the  information  recorded  is: 

(a)  the  name  of  the  current  program  -  up  to  8  characters; 

(b)  the  version  number  of  the  current  program.  From  time  to  time  program 
changes  may  be  made  sufficient  to  justify  a  new  version  number; 

(c)  the  version  number  of  the  data  file.  Each  time  UPDATE  operates,  it 
increments  this  value  by  one; 

(d)  the  current  date  and  time.  This  information  is  extracted  by  the 
computer  system  from  its  own  records.  The  method  of  storage  of  these  values  is 
seen  by  a  study  of  the  listing  of  the  subroutine  UPDATE  ; 
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(e)  the  user's  code  name.  As  with  (d) ,  this  is  provided  by  the  computer 
system; 

(f)  details  of  the  current  program  operation.  Up  to  32  characters  of 
description  may  be  supplied  by  the  user,  in  response  to  a  request  for  "details 
of  this  run". 

2.2  Data  records 

The  cartographic  data  in  a  map  file  follows  the  header  record,  and  is 
stored  in  similar  sized  data  records.  The  data  consists  of  a  sequence  of  map 
'features',  these  being  of  three  types;  corner  points,  points,  and  lines.  The 
two  types  of  point  contain  a  relatively  small  amount  of  information  and  each  is 
stored  as  one  'item',  but  lines  can  contain  an  indefinitely  large  amount  of 
information,  so  they  are  broken  down  into  a  number  of  items,  one  for  each 
coordinate  pair  of  the  line,  the  first  and  last  coordinate  pairs  being  in  rather 
different  types  of  item  from  each  other  and  from  the  'internal'  ones. 

Map  data  is  thus  broken  down  into  items,  which  may  be  from  2  to  21  real 
words  in  length,  though  the  average  item  size  is  probably  under  three  words  for 
a  typical  map.  These  items  are  stored  within  data  records  of  length  512  real 
words  (1024  integers),  and  each  type  of  item  may  carry  an  indicator  word  (denoted 
by  MK)  according  to  its  type  (see  Table  2). 

Table  2 

Meaning  of  MK  numbers 


MK 

Meaning 

0 

'Internal'  line  coordinate  pair 

1 

Last  coordinate  pair  of  a  line 

2 

Not  used 

3 

Start  of  a  line,  including  first  coordinate  pair 

4 

Map  point 

5 

Corner  point  or  'fiducial' 

6 

End  of  data  indicator 

7 

Fill  word 

Since  items  are  stored  without  being  broken,  there  is  often  one  or  more 
words  unused  at  the  end  of  a  data  record.  In  such  cases  the  first  or  only  empty 
word  has  a  'fill'  word  (MK  =>  7)  placed  in  it,  to  distinguish  it  from  a  zero. 
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The  condition  MK  =  7  is  written  and  read  by  subroutines  NSTORE  and 
NFETCH  respectively.  It  is  not  used  in  the  calling  program,  nor  need  its 
existence  be  considered  by  the  writer  of  calling  programs,  for  it  is  to  him  an 
illegal  condition.  It  is  included  in  Table  2  for  completeness  only. 

When  the  last  map  item  has  been  written,  a  MK  =  6  item  is  placed  after 
it,  to  indicate  that  the  end  of  the  data  has  been  reached. 

In  addition  to  the  coordinate  pairs,  where  appropriate,  the  item  may  also 
carry  some  information  about  the  feature,  and  this  is  done  for  star t-of-line  and 
map  point  items.  Both  such  features  may  carry  up  to  17  feature  codes,  so  the 
item  includes  NFC  ,  the  number  of  feature  codes  for  this  item,  followed  by 
NFC  real  words  containing  the  relevant  feature  codes.  Additionally  the  start-of- 
line  item  has  provision  to  record  the  number  of  steps  in  that  line  (ie  the  number 
of  coordinate  pairs  less  one).  The  number  of  steps  is  broken  down  into  two 
integer  values,  NPTS  and  NPTZ  such  that  the  total  number  of  steps  is 
NPTZ  x  10000  +  NPTS  .  In  practice,  lines  vary  in  length  from  1  to  about  50000 
steps,  so  NPTZ  seldom  exceeds  say  5.  If  the  number  of  steps  is  not  known, 
zero  should  be  inserted  in  NPTS  and  NPTZ  . 

In  the  case  of  'internal'  line  coordinate  pairs,  the  inclusion  of  the  MK 
word  would  lead  to  a  50%  increase  in  item  size,  and  since  a  map  often  consists 
mainly  of  such  items,  almost  a  50%  increase  of  data  file  size.  Thus,  for  inter¬ 
nal  points,  the  MK  word  is  not  stored.  However,  it  is  necessary  to  supply 
MK  =  0  when  internal  coordinates  are  to  be  written  by  subroutine  NSTORE  ,  and 
after  the  retrieval  of  such  data  by  NFETCH  MK  =  0  is  returned. 

The  constitution  of  the  various  types  of  item  is  shown  in  Table  3 
(section  2.2.1). 

Since  the  data  storage  method  is  almost  entirely  encompassed  by  the  sub¬ 
routine  NSTORE  (and  correspondingly  the  recovery  of  the  data  by  NFETCH)  the 
former  will  now  be  described  in  some  detail, 

2.2.1  Subroutine  NSTORE 

The  file  data  must  be  stored  in  such  a  manner  that  it  can  be  recovered 
unambiguously,  and  this  is  done  by  subroutine  NSTORE  .  This  subroutine  communi¬ 
cates  with  the  other  portions  of  the  program  via  two  COMMON  blocks,  the  first 
of  which,  COORDS  ,  must  be  suppl ied  with  all  the  information  necessary  for  the 
storage  operation,  and  the  second,  OUTCOM  ,  contains  the  output  buffer  and  its 
pointer  together  with  other  details  necessary  for  the  output  operation  of  writing  i 

to  the  storage  medium.  ! 
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The  COMMON  block  COORDS  contains  the  values: 

(a)  MK:  an  indicator  of  the  type  of  item  etc  being  stored. 

(b)  Two  real  values,  the  x  and  y  coordinates  of  the  item. 

(c)  Two  integers,  representing  the  number  of  steps  in  the  line.  These 
are  only  used  if  the  item  is  the  start  of  a  line. 

(d)  One  integer  NFC  representing  the  number  of  feature  codes  for  the 

i  tem. 

(e)  Up  to  17  real  values,  the  actual  feature  codes. 

(d)  and  (e)  are  only  supplied  if  the  item  is  a  map  point  or  the  start  of  a 
line.  The  number  of  codes  supplied  in  (e)  should  be  consistent  with  (d),  as 
NSTORE  stores  NFC  values,  irrespective  of  whether  these  have  been  deliberately 
supplied  or  not. 

The  COMMON  block  0UTC0M  contains  in  particular  a  buffer  pointer  which 

points  to  the  next  free  word  in  the  output  buffer. 

The  items  placed  into  the  buffer  by  NSTORE  are  made  up  as  shown  in 
Table  3. 

Table  3 


Types  of  item  stored 


MK 

R 

X 

Y 

NFC 

Feature  codes 

Item  size  (words) 

0 

/ 

/ 

2 

1 

/ 

/ 

/ 

3 

2 

Illegal 

3 

/ 

/ 

/ 

/ 

NFC  codes 

4  +  NFC 

4 

/ 

/ 

/ 

/ 

NFC  codes 

4  +  NFC 

5 

/ 

✓ 

/ 

3 

6 

1 

7 

</) 

(0 

The  only  legitimate  MK  values  which  may  be  supplied  to  NSTORE  are 
0,  1,  3  to  6,  and  all  others  cause  an  error  message  to  be  output.  NFETCH 
recognizes  MK  =  7  as  a  fill  word. 

Subroutine  NSTORE  starts  by  using  MK  ,  and  if  necessary  NFC  ,  and 
calculates  the  size  of  the  item,  as  shown  in  the  last  column  of  Table  3.  It 
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then  uses  the  buffer  pointer  and  buffer  size  to  determine  whether  sufficient  space 
remains  to  write  the  item  to  the  buffer.  If  insufficient  space  remains,  it 
'closes'  the  buffer  and  writes  it  to  the  output  device.  Closing  the  buffer  con¬ 
sists  of  examining  it  to  see  if  any  empty  space  remains,  and  if  so,  an  R  value 
corresponding  to  MK  =  7  is  placed  in  the  first  or  only  empty  space.  R  is  a 
real  number,  made  up  of  the  addition  of  MK  +  1048576.0  .  After  writing  the 
buffer  to  the  output  device,  the  buffer  is  cleared  and  its  pointer  reset. 

In  the  event  that  sufficient  space  remains  in  the  buffer,  or  a  buffer  has 
just  been  written  and  cleared,  the  item,  as  described  in  Table  2,  is  then  loaded 
into  the  buffer.  Much  of  this  operation  is  straightforward,  but  two  points  need 
elucidation. 

(a)  It  is  necessary  for  NSTORE  to  write  R  (MK  +  RMARK)  rather  than 
merely  MK  ,  to  avoid  ambiguity  when  the  data  is  later  recovered  by  NFETCH  . 

If  only  MK  were  used,  the  possibility  of  confusion  would  occur  when  line-points 
were  being  recovered.  At  any  time,  it  will  not  be  known  whether  the  next  item  is 
another  line-point  (in  which  case  the  next  word  will  have  the  value  x)  or  the 
end  of  the  buffer  (MK  =  7)  or  the  end  of  the  line  (MK  =  1).  Since  1.0  and  7.0 
are  legitimate  values  for  x  ,  confusion  could  arise.  Hence  MK  has  RMARK 
added  to  it,  RMARK  being  large  enough  to  exceed  any  allowed  value  of  x  yet 
small  enough  for  the  addition  of  MK  ( 1  to  7)  to  allow  subsequent  recovery  of  MK 
by  NFETCH. 

The  user  of  subroutines  NFETCH  and  NSTORE  need  not  be  concerned  with  the 
details  of  the  construction  of  R  since  it  is  handled  internally  by  these  two 
subroutines  and  no  others. 

(b)  The  values  of  the  three  integers  NFC  ,  NPTS  and  NPTZ  are  combined 
into  one  real  word.  This  is  done  for  reasons  of  economy;  as  NFC  cannot  exceed 
17,  and  NPTZ  is  in  practice  seldom  greater  than  say,  5  to  10,  these  two  small 
integers  can  conveniently  be  packed  into  one  16-bit  word,  which  in  fact  allows 
NPTZ  values  up  to  1023. 

3  DESCRIPTION  OF  PROGRAMS  FOR  CARTOGRAPHIC  DATA  FILES 
3. 1  Brief  description  of  the  programs 

Map  files  must  first  be  created;  they  may  then  be  modified  in  various  ways; 
they  can  be  transformed  into  other  forms  of  file;  and  they  can  be  examined  in 
several  ways.  Programs  have  been  written  to  perform  a  number  of  these  operations. 
This  section  describes  the  programs  briefly;  a  fuller  description  and  program 
listing  appears  later. 
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3.1.1  Programs  to  create  map  files 

(a)  MAP. FAB  This  program  allows  the  user  to  create  a  map  file  manually, 
by  keying  all  the  details  into  a  VDU  .  Its  main  use  is  for  such  purposes  as  to 
form  an  'empty'  map  -  xe  merely  four  corner  points,  or  to  create  a  grid  of  hori¬ 
zontal  and  vertical  lines,  although  one  could,  laboriously,  key  in  a  whole  map 
if  the  data  were  available  in  coordinate  form. 

(b)  MAP. IN  This  program  creates  a  map  file  from  an  external  'exchange 
file'  produced  elsewhere  and  stored  on  magnetic  tape. 

(c)  IM. CONTOUR  This  program  is  named  in  the  ’ IM. '  series,  because  it 
uses  an  image  file  as  its  starting  or  input  material,  extracts  cartographic  lines 
from  it,  and  forms  a  map  file.  The  program  is  to  be  described  in  a  separate 
document,  and  has  been  described  in  a  general  way  in  Ref  2,  so  will  not  be 
detailed  in  this  document.  However,  it  is  at  present  one  of  the  main '■sources  of 
map  files. 

3.1.2  Programs  to  alter  or  transform  map  files 

(a)  MAP. HD RED  It  is  sometimes  necessary  to  inspect,  alter  or  add  to  the 

information  contained  in  the  header  record,  xe  edit  the  header,  and  this  may  be 
done  by  use  of  MAP.HDRED  . 

(b)  MAP. TRANS  Map  files  created  by  IM. CONTOUR  are  in  the  pseudo¬ 
projection  of  image  file  units.  With  the  aid  of  a  suitable  transformation  matrix, 
MAP. TRANS  can  be  used  to  transform  such  a  map  file  to  one  in  a  known  map  pro¬ 
jection,  such  as  British  National  Grid. 

(c)  MAP. OUT  This  program  uses  a  map  file  as  input  and  creates  an 
'exchange  file'  on  magnetic  tape,  which  may  then  be  despatched  elsewhere,  or 
could  of  course  be  re-read  by  MAP. IN  . 

3.1.3  Programs  to  examine  map  files 

(a)  MAP. LIST  This  program  displays  to  the  user,  on  a  VDU  ,  some  details 
of  the  header,  and  the  data,  item  by  item,  all  in  alphanumeric  form.  To  avoid 
the  display  of  large  quantities  of  data,  the  user  may  make  some  degree  of 
selection  of  what  is  to  be  displayed. 

(b)  MAP.TEK  A  map  file  may  be  graphically  displayed  on  a  Tektronix  4014 
device  by  means  of  program  MAP.TEK  . 

(c)  MAP.HPLOT  This  program  allows  a  map  file  to  be  plotted  on  a  Hewlett- 
Packard  X-Y  plotter.  Both  MAP.HPLOT  and  MAP.TEK  allow  the  user  some  degree 
of  control  over  the  material  displayed  and  the  size  of  the  display. 
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3.1.4  List  of  programs 

In  alphabetical  order,  the  programs  above  are: 

Name  Section  Description 


IM. CONTOUR 

3.1.1c 

Ref  3 

MAP. FAB 

3.1.1a 

3.2. 1 

MAP.HDRED 

3.1.2a 

3.2.2 

MAP . HP LOT 

3.1.3c 

3.2.3 

MAP. IN 

3.1.1b 

3.2.4 

MAP. LIST 

3.1.3a 

3.2.5 

MAP. OUT 

3.1.2c 

3.2.6 

MAP.TEK 

3.1.3b 

3.2.7 

MAP. TRANS 

3.1.2b 

3.2.8 

The  above  list  of  programs  is  not  intended  to  be  exhaustive,  but  includes 
all  of  those  so  far  written.  There  may  be  a  need  for  including  additional 
features  in  these  programs  or  writing  of  new  ones  for  other  tasks.  For  example, 
it  may  soon  be  useful  to  be  able  to  create  map  files  by  digitization  from  a  map 
or  other  diagram  using  the  graphics  tablet  at  present  available. 

3.2  Detailed  description  of  the  programs 

With  the  exception  of  IM. CONTOUR  ,  which  is  to  be  described  elsewhere,  the 
programs  are  next  described  in  more  detail,  the  aim  being  to  aid  potential  users 
in  understanding  how  to  operate  them.  The  program  listings  are  provided  at  the 
end  of  this  document,  and  these  may  be  referred  to  whilst  studying  the  following. 

The  programs  make  use  of  a  number  of  subroutines  which  have  been  placed 
into  a  'library',  and  are  described  in  section  4. 

3.2.1  Program  MAP. FAB 

Program  MAP. FAB  allows  the  user  to  create  a  map  file,  by  supplying  data 
from  a  keyboard. 

From  the  user's  viewpoint,  its  operation  is  as  follows: 

*  The  program  name  and  writing  date  is  output. 

*  An  output  filename  ( ie  the  file  to  be  created)  is  requested. 

*  Permanent  header  detail  is  requested. 

*  Description  of  the  current  run  is  requested.  The  header  record  is  then 
written  to  the  storage  medium.  If  further  header  detail  is  to  be  supplied, 
this  is  done  by  subsequent  use  of  program  MAP.HDRED  . 

*  A  'MK'  number  is  requested.  Subsequent  action  depends  on  the  reply  given. 

(a)  MK  ■  6  causes  the  output  file  to  be  closed  and  the  program  terminates. 
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(b)  Illegal  MK  numbers  invoke  an  appropriate  message,  and  a  further 
request  for  a  MK  number. 

(c)  Legal  MK  numbers  are  followed  by  requests  for  the  data  needed  to 
create  an  item  appropriate  to  that  MK  number.  When  the  necessary 
answers  have  been  provided  the  item  is  written  by  NSTORE  and  a 
request  for  another  MK  number  is  made. 

It  is  thus  seen  that  this  program  creates  a  map  file,  with  the  header 
record  in  place  and  partially  completed,  and  followed  by  data  items,  each  checked 
for  internal  correctness,  and  terminated  by  a  MK  =  6  item.  There  is  however  no 
check  that  the  sequence  of  the  data  as  a  whole  is  correct:  for  example,  fiducials 
may  be  written  at  any  time,  and  a  line  may  be  ended  before  it  has  been  started. 
Checks  could  be  made  for  such  errors,  but  it  is  thought  desirable  to  allow  the 
user  to  include  deliberate  mistakes,  to  test  other  programs. 

3.2.2  Program  MAP. HD RED 

The  header  record  is  seen  (Table  1)  to  consist  of  a  number  of  values,  some 
integer,  some  real  and  some  characters,  in  the  first  132  (integer)  words.  Words 
201  to  800  contain  the  history  of  the  file. 

The  program  MAP.HDRED  allows  the  user  to  examine,  and  if  desired  change, 
any  of  the  information  contained  in  words  1  to  132.  If  any  change  is  made,  the 
history  is  updated  to  include  this  fact.  The  name  of  the  file  is  not  changed. 

The  operation  of  the  program  is  as  follows: 

*  The  program  announces  its  name  and  version  date. 

*  The  name  of  the  data  file  is  requested. 

*  The  question  'HEADER?*  is  asked.  If  answered  YES  the  project  and  descrip¬ 

tion  are  listed,  together  with  the  history  of  the  file,  as  stored  in 
header  words  201  onwards. 


Next  an  'ITEM'  is  requested.  Thirty-six  different  answers  are  acceptable 
as  shown  in  Table  4,  otherwise  the  request  will  be  made  again,  after  a 
message  stating  that  the  item  is  not  recognized.  If  the  user  replies 
HELP  ,  the  36  recognisable  options  will  be  listed  for  inspection.  If  END 
is  typed,  or  merely  a  carriage-return,  the  program  closes  down,  updating 
the  header  and  writing  it  back  into  the  data  file  if  necessary  (if  any 
change  has  been  made).  If  any  of  the  recognised  names  are  typed,  the  value 
of  that  quantity  will  be  displayed,  and  the  question  posed:  CHANGE  IT? 
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If  NO  is  replied,  no  action  is  taken,  and  the  program  asks  for  another 
ITEM  .  If  YES  ,  then  NEW  VALUE  is  requested,  and  the  value  provided 
is  written  into  the  header  buffer.  No  check  is  made  by  the  program  con¬ 
cerning  the  validity  of  the  value  provided,  so  the  user  must  take  care  to 
ensure  that  the  required  value  is  entered  correctly.  If  there  is  doubt,  or 
a  mistake  is  known  to  have  been  made,  it  is  possible  to  re-examine  and  alter 
the  altered  value. 

*  When  the  user  has  replied  END  to  the  query  'ITEM'  the  program  checks 
whether  any  change  has  been  made  (a  change  is  said  to  have  been  made  when 
the  reply  YES  is  given  to  the  query  CHANGE  IT?  even  if  the  same  value 
is  given  as  existed  before)  and  if  so,  it  updates  the  header  and  writes 
it  back  to  the  storage  medium.  The  file  is  then  closed. 

*  Finally  the  program  r  quests  MORE  FILES?  to  enable  the  user  to  edit  the 
header  of  another  file  should  he  so  wish. 

3.2.3  Program  MAP.HPLOT 

Program  MAP.HPLOT  is  used  to  plot  a  map  file  on  a  Hewlett-Packard  X-Y 
plotter.  The  interface  software  to  drive  that  plotter  has  been  obtained  from 
elsewhere  and  is  therefore  not  described  here.  The  various  subroutine  calls, 
eg  HP IN IT,  HPMOVE  ,  and  so  forth,  are  noted  in  the  program  listings  with 
comnent  to  indicate  their  functions. 

User  operation  is  described  below.  Before  running  the  program,  the  user 
must  assign  the  plotter  to  the  current  VDU  ,  by  the  instruction:  AS  AMLC  15 

*  Program  name  and  date  is  output  to  the  VDU  . 

*  The  pen  number  (1  to  4)  is  requested.  The  user  may  insert  pens  with 
coloured  ink  into  some  or  all  of  the  four  positions  available. 

*  The  cross  size  N  is  requested.  Corner  and  other  points  are  drawn  as 
crosses,  the  arm  length  of  which  will  be  N  units,  each  one-eighth  of  a 
millimetre.  Hence  the  reply  10  would  produce  crosses  of  arm  length 
li  mm,  ie  2}  mm  from  tip  to  tip. 

*  A  request  is  made  for  the  name  of  the  input  map  file. 

*  An  enquiry  is  made;  ’SPECIFY  CORNERS?'  If  the  reply  is  YES  ,  then 
requests  for  x  and  y  coordinates  of  the  south  west  and  north  east 
corners  {ie  lower  left,  upper  right)  will  follow,  and  the  plot  will  be 
adjusted  so  that  the  map,  based  on  the  corners  in  the  header  (CNR),  fits 
the  numbers  provided. 
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*  If  the  map  file  has  IPRJC1  =  99  ,  ie  it  is  in  image  file  units,  the  pro¬ 
gram  moves  to  the  enquiry  FRAME?  below.  Otherwise,  if  the  above  question 
were  answered  NO  ,  the  next  request  would  be  'SPECIFY  SCALE?'.  If  this 

is  answered  YES  the  program  will  refer  to  the  SCALE  in  the  header  and 
if  it  finds  no  value  there  (actually  SCALE  less  than  10)  a  suitable 
message  will  be  output.  Provided  that  a  suitable  SCALE  is  available,  the 
program  will  request  an  'OUTPUT  SCALE'  .  A  check  will  then  be  made  to 
ensure  that  the  map,  when  drawn  to  this  scale  will  fit  the  size  of  the 
plotter  and  if  unsuitable,  a  fresh  value  will  be  requested.  The  scale 
should  be  provided  as  a  real  number,  as  shown  in  section  2.1.2(c). 

*  If  the  answer  to  the  question  'SPECIFY  SCALE?'  had  been  NO  ,  then  the 
program  would  calculate  the  necessary  factor,  such  that  the  output  map 
would  just  fit  the  height  or  width  of  the  plotter,  whichever  is  the  limit¬ 
ing  factor. 

Thus  the  user  has  a  choice  of  three  types  of  map  size,  or  two  if  there  is 
no  value  for  SCALE  in  the  header. 

*  An  enquiry  'FRAME?'  is  made.  If  the  response  is  YES  ,  the  plotter  will 
draw  a  frame  round  the  map,  ie  a  line  joining  the  four  corners,  in  sequence 
and  returning  to  the  starting  corner. 

Following  these  questions,  the  X-Y  plotter  will  then  draw  the  map,  under 
computer  control,  using  the  data  from  the  input  file.  When  the  plot  is 
complete; 

*  the  query  'ANY  MORE  FILES?'  is  made.  If  YES  ,  a  return  is  made  to  the 
first  question;  if  NO  ,  the  program  stops. 

3.2.4  Program  MAP. IN 

For  purposes  of  data  exchange  or  transfer  between  users  in  different 
establishments,  a  data  'exchange  format'  has  been  defined,  and  this  is  described 
in  Appendix  A.  One  or  more  maps  may  be  contained  on  such  a  magnetic  tape.  The 
program  MAP. OUT  converts  data  files  into  exchange  format  data  on  magnetic  tape, 
and  MAP. IN  correspondingly  creates  one  or  more  map  format  data  files  from  an 
exchange  format  magnetic  tape. 

Before  the  program  MAP. IN  is  run,  the  user  must  first  allocate  or  'assign' 
a  magnetic  tape  drive  unit  for  his  use,  and  load  the  magnetic  tape  containing  the 
exchange  format  data  onto  that  drive  unit.  The  program  may  then  be  run: 

*  The  program  announces  its  name  and  asks  for  the  magtape  drive  unit  number. 
The  user  must  supply  the  unit  number  previously  allocated. 
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*  The  first  ten  magnetic  tape  records  are  read  by  the  computer,  interpreted 
as  alphanumeric  characters  and  displayed  on  the  VDU  .  These  ten  records 
usually  give  some  descriptive  details  of  the  map  to  follow. 

*  The  usual  filename,  permanent  header  data  and  run  description,  are 
requested.  The  user  is  able  to  provide  suitable  replies,  often  based  on 
the  ten  records  just  supplied. 

Following  these  questions,  the  computer  then  reads  the  magnetic  tape, 
converting  the  data  into  a  suitable  form  for  the  map  file  being  created  on 
the  storage  medium  (usually  a  disc),  until  the  first  (or  only)  map  has  been 
completed. 

If  no  more  files  exist  on  the  magnetic  tape  (and  this  is  recognized  by  the 
presence  of  a  MK  =  30  or  an  end-of-file)  then  the  tape  is  rewound  and  the 
program  stops.  If,  however,  another  file  (or  files)  is  present,  the  ten 
character  records  are  displayed  on  the  VDU  ,  and  the  program  continues 
from  that  location. 

Thus,  one  magnetic  tape  in  exchange  format  may  result  in  the  formation  of 
several  data  files.  Each  of  these  data  files  will  have  only  a  small  amount  of 
header  record  detail,  consisting  of  the  file  descriptive  character  words,  the 
update  information,  and  the  map  corners  (CNRs)  derived  from  the  corner  points. 

If  it  is  desired  to  include  further  header  information,  eg  the  geographical 
corners,  or  the  scale,  which  may  have  been  provided  by  the  ten  lines  of  text 
message,  then  MAP.HDRED  may  be  used  to  enter  this. 

The  exchange  format  allows  one  or  more  records  per  magtape  block.  Both 
MAP. IN  and  MAP. OUT  have  been  written  for  the  case  of  one  record  per  block, 
which  has  been  suitable  for  the  data  transfers  so  far  performed.  Should  it  be 
necessary  to  use  multi-record  blocks,  some  program  modification  would  be  necessary. 

3.2.5  Program  MAP. LIST 

The  contents  of  a  map  data  file  may  be  inspected  in  the  form  of  characters 
and  numbers  on  a  VDU  by  means  of  program  MAP. LIST  .  When  run: 

*  The  program  announces  its  name  and  date. 

*  A  map  filename  is  requested  and  the  file  opened. 

*  The  question  HEADER?  is  posed.  If  the  user  replies  NO  ,  the  program 
proceeds  to  the  data,  but  if  YES  ,  a  number  of  header  details  is  displayed 
on  the  VDU  .  Only  the  more  commonly  inspected  details  are  shown,  but  the 
less  frequently  used  details  may  be  seen  with  the  aid  of  MAP.HDRED  . 
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*  The  question  DATA?  is  asked.  If  the  reply  is  NO  ,  the  program  moves  on 
to  a  later  question,  but  if  YES  the  question  INTERIOR  LINE  POINTS? 
follows.  If  this  is  answered  YES  ,  then  all  of  the  map  data  items  are 
displayed,  in  sequence,  on  the  VDU  .  The  reply  NO  would  cause  the 
interior  line  point  items  (MK  =  0)  to  be  omitted,  which  is  often  convenient. 

*  Following  the  data,  the  program  asks  MORE  FILES?  If  YES  ,  a  return  is 
made  to  the  first  question  in  the  program,  but  if  NO  the  program  stops. 

3.2.6  Program  MAP. OUT 

Program  MAP. OUT  performs  the  reverse  operation  of  that  done  by  MAP. IN  , 
and  reference  may  usefully  be  made  to  the  description  of  the  latter  program,  in 
section  3.2.4. 

MAP. OUT  is  used  to  form  a  magnetic  tape,  in  exchange  format,  containing 
one  or  more  map  files. 

Before  MAP. OUT  is  run,  the  user  must  assign  a  magnetic  tape  drive  unit 
for  his  use,  and  load  onto  that  unit  the  tape  on  which  the  exchange  data  is  to 
be  written.  The  program  may  then  be  run. 

*  The  program  provides  its  name  (on  the  VDU)  and  requests  the  magtape  drive 
unit  number. 

*  The  name  of  the  input  map  data  file  is  requested,  and  the  named  file  then 
opened . 

*  Ten  header  lines  are  requested.  The  user  must  supply  ten  lines  of  text, 
from  0  to  120  characters  per  line.  Since  the  exchange  format  does  not 
contain  headers,  it  is  usual  to  write  out  some  of  the  header  information 
as  text. 

The  data  is  then  converted  to  exchange  format  and  written  onto  the  magtape. 
When  the  process  is  completed  the  program  asks  ANY  MORE  FILES?  If  YES  , 
the  name  of  the  next  file  is  required,  and  the  above  steps  continue.  If 
NO  ,  MK  =  30  and  end-of-file  markers  are  placed  on  the  tape,  which  is 
then  rewound,  and  the  program  stops. 

3.2.7  Program  MAP.TEK 

Program  MAP.TEK  is  used  (in  a  manner  analogous  to  MAP.HPLOT)  to  display 
a  map  file  on  a  Tektronix  4014  graphics  display  unit,  with  hard  copy  capability. 

As  with  HPLOT  ,  the  software  which  drives  the  display  has  been  provided  from 
elsewhere. 
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For  operation,  the  user  must  first  assign  the  display  to  his  use  by  making 
the  instruction: 

AS  AMLC  14  TRANHS  413 
with  corresponding  un-assignment: 

UN  AMLC  14 

after  the  program  use  has  been  completed.  The  assignment  instruction  may  be 
changed  at  some  later  date,  eg  the  AMLC  line  may  no  longer  be  number  14. 

When  the  program  is  run: 

*  it  provides  the  program  name  and  date  on  the  VDU, 

*  it  asks  if  the  user  wishes  to  'flash'  the  screen,  ie  clear  it  of  its 
present  contents, 

*  it  requests  CROSS  SIZE  ,  N  .  All  map  points  are  drawn  as  crosses,  having 
four  arms,  each  of  length  N  units.  These  units  are  related  to  the  size 
of  the  display  screen  but  are  approximately  0. 1  mm.  Hence  a  reply  of 

N  ■  10  would  provide  crosses  with  arms  of  length  about  1  mm,  ie  2  mm 
from  tip  to  tip. 

*  the  filename  is  requested  from  the  user,  and  when  supplied  the  map  file 
is  opened, 

*  SPECIFY  CORNERS?  is  asked.  If  the  reply  is  YES  the  user  is  asked  for 
the  x  and  y  coordinates  of  the  south  west  and  north  east  corners. 

Values  should  be  in  screen  units,  as  mentioned  above,  and  normally  these 
would  be  provided  in  the  range  0  to  5000.  It  is  possible  to  use  specified 
corners  to  fit  a  map  file  to  a  displayed  image  file.  If  the  user  replies 
NO  the  program  refers  to  the  corner  points  (CNR)  in  the  header  record  of 
the  file,  and  calculates  a  scaling  factor  such  that  the  map  is  drawn  as 
large  as  possible  consistent  with  (a)  scales  in  the  x  and  y  direction 
are  the  same  and  (b)  the  map  fits  entirely  into  the  screen. 

*  The  program  asks  FRAME?  If  YES  ,  four  straight  lines  are  drawn,  clock¬ 
wise  from  corner  to  corner,  creating  a  frame  round  the  map. 

The  map  file  is  then  drawn  on  the  screen,  feature  by  feature,  until  it  is 
complete.  The  file  is  then  closed. 

*  ANY  MORE  FILES?  is  requested.  YES  causes  a  return  to  the  initial 
'FLASH  THE  SCREEN?'  question,  ready  for  another  map  file  to  be  displayed, 
whilst  NO  causes  the  program  to  stop. 
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If  the  corner  (CNR)  values  in  the  header  record  of  a  file  are  either  2ero 
(perhaps  none  has  been  provided)  or  corrupt,  an  error  message  will  be  displayed 
and  that  file  closed.  In  these  circumstances  it  is  possible  to  use  program 
MAP.HDRED  to  insert  or  correct  the  CNR  values.  The  CNR  values  can  be 
selected  by  examination  of  the  fiducial  points  (the  first  four  data  features  of 
the  file)  using  program  MAP. LIST  . 

3.2.8  Program  MAP. TRANS 

Map  data  files  produced  from  image  files  by  means  of  programs  such  as 

IM. CONTOUR  are  not  in  a  cartographically  recognised  map  projection,  but  are  in 

a  coordinate  system  related  to  the  initial  image.  This  system  has  been  referred 
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to  as  Space  Oblique  Mercator  .  The  program  MAP. TRANS  may  be  used,  together 
with  a  suitable  transformation  file,  to  convert  the  map  file  to  a  coordinate 
system  such  as  British  National  Grid. 

The  construction  and  use  of  a  transformation  file  is  described  by 
J.R.  Williams^.  Briefly  the  process  is  as  follows.  Several  locations  are 
identified  on  the  image  and  the  corresponding  National  Grid  coordinates  are 
obtained  from  maps.  The  program  IM.MATRIX  is  then  used  to  convert  the  relative 
location  data  into  a  transformation  file  containing  six  numerical  values,  T1 
and  T2  referring  to  the  translation  of  the  coordinates  and  A1 I  to  A22 
referring  to  the  shear  and  rotation  of  the  scene. 

The  program  MAP. TRANS  uses  a  transformation  file  to  convert  a  map  file 
in  image  units  (IPRJC1  =  99)  to  another  map  file  in  British  National  Grid 
coordinates  (IPRJC1  =  1)  .  Each  output  coordinate-pair  X0  ,  Y0  is  calculated 
from  the  input  coordinates  XI  ,  YI 

X0  =  (XI  -  Tl)  x  All  1  +  (YI  -  T2)  x  All 2 
Y0  =  (XI  -  Tl)  x  AI2 1  +  (YI  -  T 2)  x  AI22 

where  All  1  to  AI22  is  the  inverse  matrix:- 

AIJ1  =  A22/D 

All  2  =  -A12/D 

AI21  =  -A21/D 

AI22  =  All/D 

and  D  -  (A1 1  x  A22  -  A21  x  A12)  . 

The  minus  signs  for  All 2  and  AI21  are  needed  because  an  image  file  has  its 
Y  direction  downwards,  ie  in  the  opposite  direction  to  a  map,  or  map  file. 
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It  should  be  noted  that  the  transformation  matrix  must  have  been  calculated 
from  control  points  whose  image  coordinates  refer  to  the  same  portion  of  the 
image  which  gave  rise  to  the  map  file.  It  is  not  uncommon  for  a  user  to  select 
a  subscene  of  an  image  (using  program  IM.SUB).  A  map  file  made,  by  means  of 
IM. CONTOUR  ,  from  such  a  subscene  could  not  employ  the  matrix  derived  for  the 
whole  scene.  If  this  were  attempted,  the  National  Grid  map  would  be  displaced 
in  northing  and  easting. 

Operation  of  MAP. TRANS  is  as  follows: 

*  The  program  announces  its  name  and  date. 

*  A  matrix  filename  is  requested.  This  file  is  then  opened  by  the  computer, 
relevant  values  read  and  the  file  closed. 

*  The  input  map  filename  is  requested  and  the  file  opened.  A  check  is  made 
that  the  file  is  in  the  required  pseudo-projection.  If  it  is  not,  the  file 
is  closed  and  a  return  is  made  for  a  new  matrix  filename.  Assuming  a 
suitable  file  has  been  selected,  an  output  filename  is  requested,  and 
opened;  an  update  is  made,  including  a  request  for  a  description  of  the 
run;  a  number  of  calculations  is  made  and  the  geographical  corner  points 
(the  corners  expressed  as  kilometres  on  the  National  Grid  projection)  are 
printed  out,  and  placed  into  a  header  buffer,  together  with  other  header 
data.  The  output  header  record  is  then  written  to  the  output  medium,  and 
the  data  of  the  file  is  converted  to  the  output  projection,  and  stored. 

When  all  data  has  been  transformed,  the  input  and  output  files  are  closed. 

*  MORE  FILES?  is  asked.  If  YES  the  program  returns  to  the  request  for 
a  matrix  filename;  if  NO  ,  it  stops. 

4  SUBROUTINE  LIBRARY 

All  of  the  programs  described  in  section  3  call  subroutines.  Inspection 
of  the  program  listings  reveals  that  some  of  the  subroutines  are  called  by  more 
than  one  program.  In  such  cases  it  is  convenient  to  form  a  'library'  of  sub¬ 
routines,  which  are  accessible  without  being  re-listed  and  re-compiled  for  each 
program. 

The  subroutines  used  by  the  map  file  programs  and  placed  into  the  library 
include  (in  alphabetical  order): 


•  V  X 
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Subroutine  name  Description  in  section  number 


NCLOSI 

4.2.10 

NCLOSO 

4.2.6 

NFETCH 

4.2.9 

NHDIN 

4.2.8 

NOPENI 

4.2.7 

NOPENO 

4.2.1 

NSTART 

4.2.4 

NSTORE 

4.2.5 

PERMHD 

4.2.2 

UPDATE 

4.2.3  . 

Several  of  these  subroutines  are  concerned  with  the  input  and  output  of 
data,  which  is  often  done  by  means  of  routines  which  are  not  standard  FORTRAN 
ones.  In  the  subroutine  listings,  the  comment  includes  references  to  PDR  3106 
and  PDR  3110  together  with  page  numbers,  on  some  occasions  when  non-standard 
subroutine  calls  are  made.  These  PDR  numbers  refer  to  the  host  computer  manu¬ 
facturer's  handbooks  and  are  detailed  in  numbers  6  and  7  respectively  in  the 
reference  list  of  this  document.  In  passing,  it  may  be  noted  that  if  it  were 
necessary  to  modify  the  programs  for  use  on  a  different  host  computer,  many  of 
the  alterations  would  occur  in  these  subroutines  at  these  input  and  output  calls. 

4. 1  Subroutine  groups 
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The  subroutines  listed  above  belong  to  one  or  other  of  two  groups;  for 
'output'  or  'input'; 

(a)  subroutines  used  when  a  map  file  is  being  written  or  'output': 

NOPENO,  UPDATE,  NSTART,  NSTORE,  NCLOSO  .  Additionally,  if  the  map  file  is  being 
created  from  a  source  which  is  not  itself  a  map  file,  PERMHD  may  be  used, 
usually  before  UPDATE  ; 

(b)  subroutines  used  when  a  map  file  is  being  read  or  'input':  NOPENI, 
NHDIN,  NFETCH,  NCLOSI  . 

The  subroutines  within  each  group  are  called  in  the  sequence  given  above. 
If  both  input  and  output  occurs  (as  for  instance  in  program  MAP. TRANS)  then 
calls  from  the  two  groups  may  be  intermingled,  but  the  sequence  within  each 
group  will  still  be  preserved.  Programs,  such  as  MAP.HDRED  which  edit, 
ie  change,  an  existing  map  file,  may  present  certain  problems.  Since  only  one 
such  program  exists  at  present,  its  special  subroutines,  for  opening  and  closing 
files  for  both  input  and  output,  and  for  writing  a  header  record  back  into  an 
existing  file,  are  not  included  in  the  library. 


As  far  as  possible  the  subroutines  contained  within  group  (a)  and  group  (b) 
have  been  designed  to  form  self-contained  sets,  needing  the  minimum  of  external 
information.  Subroutine  parameters  are  not  used,  and  the  transfer  of  information 
is  by  means  of  COMMON  blocks.  For  group  (a),  the  output  routines,  the  COMMON 
blocks  are  OUTCOM,  COORDS,  HEADER  and  FILENC  ,  The  input  routines  use  COMMON 
blocks  INCOM,  COORDS  and  HEADER. 

The  method  of  use  of  the  various  subroutines  can  be  well  seen  in  such 
programs  as  MAP. FAB  which  uses  all  of  group  (a)  in  a  simple  manner,  and 
MAP. LIST  which  uses  all  of  the  routines  of  group  (b) . 

4.2  Subroutine  descriptions 

The  subroutines  are  described  in  detail  in  Appendix  B,  being  placed  in  the 
sequence  indicated  by  the  groupings  of  section  4.1.  These  descriptions  may 
advantageously  be  studied  together  with  the  library  computer  listings  at  the  end 
of  the  Report. 

The  following  facts  are  common  to  all  of  the  subroutines  and  are  therefore 
not  included  in  Appendix  B  for  each: 

all  are  written  in  Prime  FORTRAN 

all  were  written  by  the  author  of  this  paper 

none  has  any  parameters. 
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Appendix  A 
EXCHANGE  FORMAT 

A  cartographic  data  exchange  format  has  been  defined  by  the  Experimental 
Cartography  Unit  of  the  Natural  Environment  Research  Council.  Since  no  freely 
available  description  exists,  this  Appendix  gives  full  details  of  the  format. 

The  purpose  of  the  exchange  format  is  to  allow  transfer  of  cartographic 
data  between  different  establishments,  and  in  particular  between  different  types 
of  computer.  The  method  selected  as  most  suitable  to  cover  the  widest  range  of 
equipment  is  to  write  characters  (not  numbers)  onto  magnetic  tape. 

The  format  description  is: 

*  Data  to  be  written  on  magnetic  tape:  7  or  9  track:  any  of  the  standard 
densities;  parity  even  or  odd;  any  character  set.  This  allows  a  wide 
choice,  but  the  options  used  must  be  defined  on  a  separate  document. 

*  The  tape  may  contain  one  or  more  sets  of  data,  'maps'. 

*  Each  map  consists  of  a  number  of  records,  each  containing  120  characters. 

*  Each  block  written  tr-  the  magnetic  tape  must  contain  the  same  (integer) 
number  of  records.  Frequently  one  record  per  block  is  used. 

*  Each  map  consists  of  the  sequence: 

ten  records  containing  text 
an  unspecified  number  of  data  records 

one  MK  =  6  record  (see  below)  to  indicate  the  end  of  the  map. 

*  After  the  final  map  has  been  written,  including  its  MK  =  6  record,  a 
MK  =  30  record  is  added  to  indicate  the  end  of  the  data. 

*  The  text  records  may  contain  any  message  or  other  information,  in  character 
form.  Unused  portions  of  records  should  be  filled  with  space  characters  or 
nulls. 

*  Data  records  are  of  four  possible  types: 

MK  =  3  record:  start  of  line 
MK  =  4  record:  map  point 
MK  =  5  record:  corner  point 

record  containing  line  points  after  the  start  of  the  line. 
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Appendix  A 


Table  A!  lists  the  contents  of  the  MK  records. 

Table  A1 

Contents  of  MK  records 


Notes: 

(1)  Two  zero  or  space  characters. 

(2)  Six  zero  or  space  characters. 

(3)  For  MK  =  5  NFC  would  normally  be  zero. 

(4)  Up  to  17  six-character  feature-codes,  followed  by  space  or  zero 
characters.  Codes  in  the  range  -99999  to  99999. 

(5)  This  portion  of  the  record  to  be  filled  with  zero  or  space  characters. 

(6)  For  all  integers,  leading  zeros  may  be  replaced  by  spaces. 

(7)  All  numbers  are  integers. 

(8)  x  and  y  coordinates  may  be  expressed  in  any  units  but  should 
preferably  be  defined  in  an  accompanying  document. 

If  MK  =  5  records  are  supplied  (corner  points)  there  must  be  four  of  them, 
representing  respectively  the  SW,  NW,  NE  and  SE  corners  of  the  map,  and 
they  must  appear  before  all  other  data  records,  but  after  the  ten  lines  of 
text. 

For  lines,  after  the  first  point  has  been  defined  by  a  MK  *  3  record,  the 
remaining  points  are  stored  as  changes  of  x  and  y  from  the  preceding 
value.  The  change  must  never  exceed  the  range  -49  to  +49  units,  and 
should  a  larger  change  occur,  extra  points  must  be  interpolated.  Each 
change  then  has  the  value  50  added  to  it,  to  bring  it  into  the  range 
1  to  99,  the  value  zero  being  reserved  to  indicate  that  the  end  of  the 
line  has  been  reached. 


037 


Appendix  A 


27 


Line  point  records  consist  of  30  sets  of  four  characters,  two  numeric 
characters  for  the  x  change  and  two  for  the  y  change.  After  the  last  point 
of  the  line,  there  must  be  a  set  of  four  zeros  to  indicate  the  end  of  the  line, 
even  if  this  necessitates  a  new  record.  The  remaining  unused  portion  of  a  record 
may  contain  any  character. 
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B.  1  Subroutine  NOPENO 
Summary 

Use  of  COMMON 


Local  storage  used 


Appendix  B 

SUBROUTINE  SPECIFICATIONS 

-  NOPENO  is  used  to  open  a  new  map  file  for  output. 

-  uses  /OUTCOM/.  Output  arguments: 

IOTDEV  output  device  line  or  unit  number 
IOTSZ  size  of  output  buffer  in  integer  words. 

-  21  integers;  1  logical  variable;  11  bytes  for  characters; 

Prime  AgKEYS  list. 

Subordinate  subprograms  -  Prime  Applications  Library  routine  OPVPgA 
Operation  -  the  COMMON  values  IOTDEV  =  6  and  IOTSZ  =  1024 

are  provided  for  use  by  following  subroutines.  The  applications  library  OPVPjJA 
is  called,  operating  in  the  mode  A$0VAP  ,  the  effect  of  which  is  as  follows.  The 
user  is  asked  for  a  filename  which  may  be  a  Prime  'tree  name'.  OPVP^A  attempts  to 
open  a  new  file  of  this  name,  for  output.  If  the  opening  operation  is  not  success¬ 
ful,  for  example  an  illegal  tree  name  is  provided,  a  fresh  name  is  requested,  and 
this  is  repeated  until  a  successful  operation  is  achieved.  If  the  name  of  an 
existing  file  is  provided,  the  user  is  asked  if  he  wishes  to  alter  it.  If  NO  , 
a  new  name  is  requested,  but  if  YES  the  user  is  asked  whether  he  wishes  to  over¬ 
write  the  existing  file  or  to  extend  it.  The  correct  usage  for  map  files  is  either 
to  choose  a  new  filename  or  else  to  overwrite  an  existing  but  unwanted  file. 

B.2  Subroutine  PERMHD 

-  the  subroutine  PERMHD  clears  the  header  buffer, 
obtains  the  permanent  header  details  from  the  user  and 
places  them  in  the  header  buffer. 

-  uses  /HEADER/  and  /OUTCOM/ 

Input  argument  IOTSZ  in  /OUTCOM/  is  used. 

All  words  in  /HEADER/  have  values  output. 

-  1  integer,  30  characters  in  FORMAT  statements. 

-  none 

-  The  value  IOTSZ  is  used  to  determine  the  number  of 


Summary 


Use  of  COMMON 


Local  storage  used 
Subordinate  subroutines 


Operation 

words  of  the  COMMON  block  /HEADER/  ,  and  these  are  all  cleared  to  zero.  The 
user  is  then  asked  to  supply  the  permanent  header  details,  which  consist  of  the 
project  name,  IPROJ  ,  up  to  12  characters  and  the  project  description,  IDE SCR  , 
up  to  32  characters. 

PERMHD  is  usually  only  called  when  a  map  file  is  being  created.  Subsequent 
alterations  to  the  map  file  do  not  necessitate  a  call  of  PERMD  . 

B.3  Subroutine  UPDATE 


Sumnary 


-  UPDATE  is  used  to  maintain  the  file  log  in  the  header  buffer. 
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Use  of  COMMON 


Local  storage  used 
Subordinate  subprograms 

Operation 


-  uses  /FILENC/  and  /HEADER/ 

Input  arguments:  /FILENC/ 

NPR0G(4)  name  of  the  program  calling  UPDATE 
NVERS  version  number  of  that  program 
Output  arguments:  modifications  to  /HEADER/  words 
20!  to  800. 

-  17  integers,  25  characters  in  FORMAT  statements 

-  Prime  file  system  subroutine  TIMDAT 

-  Map  files  contain  a  log  in  the  header. 


Each  log  entry  is  40  words  long,  and  there  is  space  for  15  log  entries  in  header 
words  201  to  800.  Each  time  UPDATE  is  called  it  moves  the  existing  14  most 
recent  entries  (words  201  to  760)  into  the  positions  for  entries  2  to  15 
(words  241  to  800)  and  the  15th  entry,  if  it  existed,  is  overwritten  and  hence 
lost.  A  new  entry  is  then  provided  in  heauer  words  201  to  240  as  follows.  The 

name  and  version  number  of  the  calling  program  (which  must  be  placed  in 

COMMON /frlLENC/  by  that  program)  are  moved  into  header  words  201  to  205.  The  log 
sequence  number,  in  header  word  206  is  obtained  from  the  sequence  number  of  the 

previous  log,  in  word  246,  incremented  by  one.  For  the  first  UPDATE  call 

word  246  is  zero  so  the  first  log  number  is  made  1  ,  as  required.  The  Prime  sub¬ 
routine  TIMDAT  is  then  called,  and  this  moves  information  about  the  date,  time 
and  user-code  into  the  dimension  IARRAY  ,  from  which  it  is  moved  into  the  appro¬ 
priate  header  words.  A  request  is  then  made  to  the  user  for  details  of  the 
computer  run  being  performed  and  up  to  32  characters  can  be  stored. 


B.4  Subroutine  NSTART 

Summary  -  NSTART  is  used  to  write  the  header  record  of  a  map 

file,  and  prepare  the  output  buffer  and  pointer  for 
data. 

Use  of  COMMON  -  Uses  /HEADER/  and  /0UTC0M/ 

Input  arguments  -  all  of  /HEADER/  ;  IOTDEV  and 
I0TSZ  (set  by  subroutine  NOPENO) 

Output  arguments  -  IOTBUF  and  I0TPTR  ,  the  output 
buffer  and  its  pointer 

Local  storage  used  -  1  integer,  25  characters  in  FORMAT  statements 

Subordinate  subprograms  -  Prime  file  system  subroutine  PRWF$$ 
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Operation  -  The  entire  contents  of  the  COMMON  block  /HEADER/  are  transferred 
to  the  output  buffer  and  written  to  the  output  file,  from  its  start,  by  Prime 
subroutine  PRWFgg  .  The  output  file  is  on  device  line  IOTDEV  and  the  buffer 
size  is  IOTSZ  .  If  a  system  failure  occurs  whilst  the  write  operation  is  taking 
place,  an  error  message  ' PRWF$$  FAILURE  IN  NSTART'  is  sent  to  the  terminal  and 
the  program  stopped.  All  values  in  the  output  buffer  are  then  cleared  (set  to 
zero)  and  the  output  buffer  pointer,  IOTPTR  is  set  to  1  ,  ready  for  data 
to  be  placed. 

B.5  Subroutine  N STORE 

Summary  -  NSTORE  is  used  to  place  one  data  item  to  the  output  buffer, 

and  when  the  latter  is  full  to  write  it  to  the  map  file. 

Use  of  COMMON  -  Uses  /COORDS/  and  /0UTC0M/ 

Input  arguments  (a)  in  /COORDS/ 

MK  item  type  indicator 
X,Y  coordinate  pair 

NPTS,  NPTZ  number  of  'points'  in  the  item 
NFC  number  of  feature-codes  in  the  item 
FC0DE(17)  up  to  17  feature-codes 
(b)  in  /0UTC0M/ 

IOTDEV  output  file  line  number  -  set  by  NOPENO 
IOTSZ  output  buffer  size  -  set  by  NOPENO 
Output  arguments  in  /0UTC0M/ 

IOTPTR  output  buffer  pointer 
IOTBUF  output  buffer 

Local  storage  used  -  5  integers,  2  real  variables,  38  characters  in  FORMAT 

statements,  Prime  list  KEYS.F 

Subordinate  subprograms  -  Prime  file  system  subroutine  PRWFgj5 

Operation  -  Map  data  is  stored  in  the  form  of  items  (see  the  body  of  this 

document  for  details  of  items),  each  of  which  consists  of  a  MK  number  in  the 
range  0-6,  followed  by  an  appropriate  amount  of  map  data.  In  the  case  of  MK=0, 
the  MK  number  is  not  itself  stored,  but  its  presence  can  be  inferred.  Items 
are  written  to  the  output  buffer  IOTBUF  ,  which  is  written  to  the  output  storage 
medium  when  it  has  insufficient  space  for  the  next  item,  or  following  a  MK=6  . 

The  detailed  operation  of  NSTORE  is  as  follows: 

The  variable  ISWTCH  is  set  to  1  ,  to  control  the  normal  routing  of  the  sub¬ 
routine  flow.  Only  for  the  case  MK *  6  can  ISWTCH  have  any  other  value.  The 
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value  of  MK  is  checked,  and  the  illegal  values,  ie  other  than  0,  1,  3-6, 
cause  the  subroutine  to  output  the  error  message  MK  =  XX  IN  NSTORE  and  return. 
The  item  size  in  integers  ITEMSZ  is  calculated,  from  the  value  of  MK  and 
where  necessary  NFC  .  A  check  is  then  made  of  whether  sufficient  space  remains 
in  the  buffer  to  place  the  item.  If  insufficient  space  exists,  the  buffer  must 
be  written  to  the  output  storage  device,  using  the  Prime  file  systems  routine 
PRWFjSg  to  add  the  data  to  the  end  of  the  existing  data.  Before  the  buffer  is 
written,  it  is  checked  for  any  empty  words,  and  if  one  or  more  empty  spaces 
remain,  the  first  of  these  has  a  fill-word  placed  into  it,  appropriate  to  MK=7. 
This  value  of  MK  is  internal  to  the  subroutines  NFETCH  and  NSTORE  and  need 
not  be  considered  by  the  user  of  the  subroutine  -  indeed  it  is  illegal  external 
to  the  routine.  After  the  buffer  has  been  written,  all  of  its  values  are  set  to 
zero,  and  the  pointer  IOTPTR  is  set  to  1. 

If  sufficient  space  exists  in  the  buffer,  the  next  item  is  placed  into  it. 
MK  is  not  stored  as  such,  because  it  could  be  confused  with  data,  but  has  RMARK 
(1048576.0)  added  to  it.  This  value  has  been  chosen  as  it  is  large  enough  to 
avoid  confusion  with  any  possible  data  value,  yet  is  small  enough  to  allow  the 
MK  to  be  extracted  from  it.  The  other  words  in  the  item  are  placed  as  required 
by  the  MK  value,  and  may  include  some  or  all  of:  X  Y  NPTS  NPTZ  NFC  and 
up  to  17  feature  codes.  To  aid  in  placing  the  real  words,  the  buffer  is  equiva- 
lenced  to  a  real  array  OUTBUF  and  the  real  buffer  pointer  IPTR  is  used  to 
point  to  the  real  word  locations. 

For  the  case  MK  =  6  ,  which  indicates  that  the  end  of  the  map  file  has 
been  reached,  the  item  consists  only  of  one  word,  RMARK  +  MK  ,  and  this  word 
must  be  placed  into  the  output  buffer  and  written  at  once  to  the  storage  device. 
If  there  is  space  in  the  output  buffer,  the  data  can  be  placed  and  written,  but 
if  the  buffer  is  full,  the  buffer  must  first  be  written,  then  the  MK  =  6  data 
placed  and  then  a  second  'write'  take  place,  ie  a  second  call  of  PRWFJ55S  .  In 
the  case  MK  =  6  ,  ISWTCH  is  set  to  2,  and  this  allows  the  correct  routing  after 
the  PRWFjSg  call,  ie  a  return  is  not  made  following  a  PRWFgg  call,  when  a 
MK  =  6  has  been  written. 

B.6  Subroutine  NCLOSO 

Summary  -  NCLOSO  is  used  to  close  an  output  map  file 

Use  of  COMMON  -  uses  /0UTC0M/  .  Input  argument: 

IOTDEV  -  output  file  line  number 

Local  storage  used  -  1  logical  variable,  25  bytes  for  characters, 

Prime  list  A?KEYS 
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Subordinate  subprograms  -  Prime  applications  library  subroutine  CLOS$A 

Operation  -  The  Prime  subroutine  CLOSgA  is  used  to  close  the 

output  map  file  on  line  number  IOTDEV  .  If  the  CLOS$A  operation  is  unsuccess¬ 
ful  the  message  CLOSgA  FAILURE  IN  NCLOSO  is  output  to  the  terminal  and  the 
program  is  stopped. 

B.7  Subroutine  NOPENI 

Summary  -  NOPENI  is  used  to  open  a  map  file  for  input 

Use  of  COMMON  -  uses  /INCOM/  .  Output  arguments: 

INDEV  input  device  line  number 

INSZ  size  of  input  buffer  in  integer  words 

Local  storage  used  -  21  integers,  1  logical,  10  bytes  for  characters. 

Prime  AgKEYS  list 

Subordinate  subprograms  -  Prime  Applications  Library  routine  OPNPJSA 

Operation  -  the  COMMON  values  INDEV  =  5  and  INSZ  =  1024  are 

provided,  for  use  by  following  input  file  subroutines.  The  Prime  Applications 
Library  routine  OPNPgA  is  called.  This  outputs  the  request  for  an  input  file¬ 
name  at  the  terminal,  accepts  a  Prime  tree  name,  checks  it  for  validity  and  then 
attempts  to  open  an  existing  file  of  that  name.  If  the  opening  operation  is  not 
successful,  the  process  is  repeated  until  a  satisfactory  open  occurs. 

B.8  Subroutine  NHDIN 

Summary  -  reads  one  data  block  from  the  input  map  file,  and 

places  it  into  /HEADER/ 

Use  of  COMMON  -  uses  /HEADER/  and  /INCOM/ 

Input  arguments:  /INCOM/  words: 

INDEV  input  file  line  number 
INSZ  size  of  input  buffer 
Output  arguments:  /INCOM/  word: 

INPTR  -  input  buffer  pointer 
All  of  /HEADER/ 

Temporary  use:  /INCOM/ 

INBUF  input  buffer 

Local  storage  used  -  3  integers,  24  bytes  for  characters. 

Prime  list  KEYS.F 

Subordinate  subprograms  -  Prime  file  system  subroutine  PRWF$$ 


►  ♦ 
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Operation  -  The  Prime  subroutine  PRWF$$  is  called  to  read  the  next  INSZ 
words  from  the  input  map  file  into  the  input  buffer.  INSZ  has  been  set  to  1024 
by  NOPENI  ,  and  is  equal  to  one  block  of  data  on  the  Prime  disc  storage  on  which 
map  files  are  held.  NHDIN  is  intended  to  be  called  before  other  subroutines 
which  read  data  from  a  map  file,  and  it  should  therefore  read  the  first  1024  words 
of  the  file,  which  is  the  map  file  header.  If  the  PRWF$$  operation  is  unsucess- 
ful  (ICODE  not  zero)  or  the  number  of  words  read  IRNW  is  not  equal  to  INSZ  , 
NHDIN  sends  the  message  PRWFS58  FAILURE  IN  NHDIN  to  the  terminal,  and  then  stops 
program  operation. 

Following  a  successful  PRWFgg  read  operation,  the  data  in  the  input  buffer 
is  transferred  into  COMMON  /HEADER/  .  The  input  buffer  pointer  INPTR  is  set 
to  INSZ  +  1  to  ensure  that  the  following  NFETCH  call  will  be  forced  to  start 
with  a  read. 

B.9  Subroutine  NFETCH 

Summary  -  NFETCH  is  used  to  retrieve  one  map  'item'  from  the 

input  buffer,  reading  more  data  from  the  input  file  to 
the  buffer  when  the  latter  needs  replenishment. 

Use  of  COMMON  -  Uses  /COORDS/  and  /INCOM/ 

Input  arguments,  in  /INCOM/: 

INDEV  input  file  line  number 

INSZ  size  of  input  buffer 
INPTR  input  buffer  pointer 
Output  arguments: 

Some  or  all  of  /COORDS/ 

Local  storage  used  -  4  integers,  3  real  words,  37  bytes  for  characters 

Prime  list  KEYS.F 

Subordinate  subprograms  -  Prime  File  System  subroutine  PRWF$$ 

Operation  -  Map  data  is  stored  in  the  form  of  'items',  the  details 

of  which  are  controlled  by  a  MX  word.  The  subroutine  NFETCH  ,  when  called, 
retrieves  the  next  item  from  the  input  map  file,  and  places  it  into  the  relevant 
words  of  the  COMMON  block  /COORDS/  .  Items  are  read  from  the  map  file  in 
blocks  of  1024  words,  each  block  containing  many  items. 

NFETCH  starts  by  setting  all  values  of  /COORDS/  to  zero,  so  that  only 
newly  retrieved  values  are  available,  not  values  left  over  from  a  previous  call. 
If  the  input  pointer  has  a  value  less  than  the  buffer  size,  there  must  still  be 
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data  remaining  within  the  buffer,  so  this  data  is  examined.  The  next  data  word 
is  studied  to  see  if  it  contains  a  valid  MK  word,  and  if  so,  the  action 
appropriate  to  that  MK  word  is  taken.  If,  however,  the  next  word  does  not 
correspond  to  a  valid  MK  word,  the  data  must  be  an  X  coordinate,  implying 
a  MK  =  0  .  A  MK  =7  indicates  tha  Lnere  is  no  more  data  in  the  buffer,  and 
the  subroutine  path  is  then  the  same  as  for  INPTR  greater  than  INSZ  ,  ie  a 
new  buffer  of  data  is  read  from  the  data  file  by  means  of  a  PRWF$jS  call.  The 
PRWFgg  call  is  followed  by  a  check  for  successful  operation,  and  if  unsuccessful 
the  message  PRWFjS#  FAILURE  IN  NFETCH  is  sent  to  the  terminal  and  the  program 
stopped.  Following  a  successful  PRWF$j$  call,  the  input  buffer  contains  new 
data  items. 

If  the  existing  data  buffer  contains  more  data  or  if  a  new  buffer  of  data 
has  been  read  from  the  input  file,  then  the  next  item  can  be  retrieved,  the  MK 
number  being  calculated  by  subtraction  of  RMARK  ,  or  the  value  MK  =  0  obtained 
as  explained  above.  The  MK  value  indicates  which  words  are  available  and  these 
are  extracted  from  the  input  buffer,  and  placed  into  the  /COORDS/  words:  X, 

Y,  NPTS,  NPTZ ,  NFC  and  the  feature  codes  FC0DE(17).  The  input  buffer  pointer 
INPTR  and  the  real-word  buffer  pointer  IPTR  are  used  to  aid  this  process, 
pointing  to  the  relevant  locations  in  the  equivalent  integer  or  real  buffers 
INBUF  and  RINBUF  .  A  detailed  description  of  the  various  types  of  item  is  given 
elsewhere  in  this  document. 


If  the  illegal  MK  =  2  is  fourd,  the  message  MK  =  2  IN  NFETCH  is 
sent  to  the  terminal,  and  a  return  made. 


B.10  Subroutine  NCLOSI 


Summary 

Use  of  COMMON 


Local  storage  used 


-  NCLOSI  is  used  to  close  an  input  map  file 

-  Uses  /INCOM/.  Input  argument: 

INDEV  input  file  line  number 

-  1  logical  variable,  24  bytes  for  characters. 
Prime  list  A$KEYS 


Subordinate  subprograms  -  Prime  Applications  Library  subroutine  CLOS$A 

Operation  -  The  Prime  subroutine  CLOSfcA  is  used  to  close  the 

input  map  file  on  line  number  INDEV  .  If  the  CLOS$A  operation  is  unsuccessful 
the  message  CLOS$A  FAILURE  IN  NCLOSI  is  output  to  the  terminal  and  the  program 
is  stopped. 
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Table  1 


HEADER  RECORD 


Word  No. 

M  FORTRAN  format 

Name 

or  type 

1-6 

IPR0J 

6A2 

7-22 

IDESCR 

1 6A2 

23-40 

41 

IPRJC1 

Integer 

42 

IPRJC2 

Integer 

43-44 

SCALE 

Real 

45-46 

XSCALE 

Realj 

47-48 

YSCALE 

Real] 

49-50 

SCLLNG 

Real] 

► 

51-52 

SCLLAT 

Real 

53-54 

RADEA 

Real 

55-56 

RADEB 

Real 

57-80 

81 

NLIN 

Integer 

82 

NPT 

Integer 

83 

NSTEPS 

Integer! 

84 

NSTEPZ 

Integerj 

85 

MNFC 

Integer 

86 

NDIFFC 

Integer 

87-100 

101-1  16 

CNR 

8  real 

117-132 

GCNR 

8  real 

133-200 

201-800 

- 

See  below 

< 

'5  integers 

1  integer 

3  integers 

I  integer 
3A2 

1 6A2 

I I  integers 

801-1024 


ALLOCATION 

Meaning 

General  name  of  the  project 
Description  of  the  project 
Spare 

Map  projection  code 
Map  projection  subcode 
Map  scale 

X  and  Y  scale  adjustments 

Longitude  and  latitude  (degrees)  for 
specified  purposes 

Earth's  major  semi-axis  (metres) 

Earth's  minor  semi-axis  (metres) 

Spare 

Number  of  lines  in  the  file 
Number  of  points  in  the  file 

Number  of  line  points  in  the  file 

Maximum  number  of  codes  per  item 
Number  of  different  codes  in  the  file 
Spare 

Corners  of  the  map  (mm) 

Corresponding  corners,  on  the  Earth 
Spare 

15  sets  of  40  words,  to  'UPDATE' 

Program  name  and  version 
Data  file  version  number 
Date 
Time 

User's  code  name 

Details  of  the  current  process 

Spare 

Spare 


i 
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RECOGNISED  WORDS  FOR  PROGRAM  MAP . HDRED 


Keyed  word 
END 

(Carriage  return) 
HELP 


IPROJ 

IDESCR 

IPRJC 1 

IPRJC2 

SCALE 

XSCALE 

Y SCALE 

SCLLNG 

SCLLAT 

RADEA 

RADEB 

NLIN 

NPT 

NSTEPS 

NSTEPZ 

MNFC 

NDIFFC 


CNR(A.B) 

GCN(A,B) 


Effect 

The  data  file  is  closed 

The  list  of  keyed  words  is  displayed 
on  the  terminal 


\ 


\  Controls  header  word(s)  of  these  names 


) 

.  « 

Controls  header  dimension  CNR(A,B) 
Controls  header  dimension  GCNR(A,B) 
A  =  1  or  2  B  =  1  to  4  allowed. 
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Author 

NASA 

A.H.  Benny 


A.H.  Benny 

A.P.  Colvocoresses 

J.M.  Williams 

Prime  Computer  Inc. 
Prime  Computer  Inc. 


REFERENCES 

Title,  etc 

Landsat  Data  Users'  Handbook,  Appendix  C  (1975) 

Coastal  definition  using  Landsat  data  and  potential 
bathymetric  applications. 

Sixth  Annual  Conference  of  the  Remote  Sensing 
Society  (1979) 

RAE  Technical  Report  in  preparation 
Space  Oblique  Mercator. 

Photo gramme  trie  Engineering,  921-926  (1974) 

Geometric  correction  of  satellite  imagery. 

RAE  Technical  Report  79121  (1979) 

Preliminary  documentation  release  PDR  3106 

Preliminary  documentation  release  PDR  3110 


REPORTS  QUOTED  ARE  NOT  NECESSARILY 
AVAR  Ac  *  . . -vt  TMp  PUBLIC 

pr 


>  •  c. 


o  o  o 
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PROGRAMS 


C  PROGRAM  MAP. FAB 


C  PROCRAM  MAP. FAB 

C  TO  CREATE  A  MAP  FILE  BY  KEYING  IN  FROM  A  TERMINAL 
DIMENSION  IC0C42) 

LOGICAL  Q 

C0MN0N/C00RDS/MK,X,Y,NPTS,NPT2,NFC> FCODE<  1?  ) 
COMMON/F ILENC/NPROG( 4  ), NVERS 
COMMON/HEADER/I HDUM<  1  024  ) 

COMMON/OUT  C,0H/IOTDUN(  102  7  > 

•  INSERT  SYSCOM  >A*KEYS 

EQUIVALENCE  (  MK , 1C0( 1  )> 

DATA  MPROG, NVERS/8HMAP. FAB  ,  1/ 


URITE(  1  <  20B0  ) 

2000  FORMATCNAP  FAB  9-JAN-80') 

100  CALL  NOPENO 
CALL  PERMHD 
CALL  UPDATE 
CALL  NST ART 
1  DO  3  J=>1 <42 

3  ICO(J)=»0  /*CLE ARS  /COORDS/ 

URITE( 1. 1000) 

1000  FORMA  T( ' MK  i  '  > 

READ!  1-1001  )MK 
10B1  F ORMAT (16) 

IF  (MK  GT.6  OR.  MK.LT.0  OR.  (IK  ES  2  )  GOTO  990 
IF  (MK  EQ. 6  )  GOTO  800 
URITE( 1. 1002) 

1002  F ORMAT (  '  X  s  '  ) 

READ<  1 .1003  )X 

1003  FORMA  T ( F 1 0  2) 

URITEC  1  ,  1004  ) 

1004  FORMATt ’ Yi '  ) 

READ<  1 ,1003  )Y 

IF  (MK  EQ  0  OR . MK  EQ  1  . OR . MK  EQ  S)  GOTO  600 
0  UR  I TE (  1, 1003) 

1003  FORMA T ( ' NFC : '  ) 

READ( 1,1001 )NFC 
IF  (NFC  EQ  0  >  GOTO  20 
DO  10  J* 1 , NFC 
URI T£( 1, 1008) 

1008  FORMAT ( ' FC  :  '  > 

10  READ<  1 , 1003  )FCODE< J  ) 

20  IF  (MK  EM)  GOTO  8B0 
URI TE( 1, 1806  ) 

1006  F ORMA T ( ' HP T S i '  ) 

READ  (  1,  1001  )NPTS 
URITEC  1, 1007) 

1007  FORMAT( ' NPT Z: ‘ ) 

READ( 1,1001 )NPTZ 

800  CALL  NSTORE 

IF  (MK  HE  6 )  GOTO  1 
CAt  L  NCLOSO 

Q»  SNO*A( ' MORE  FILES' , IB, ASNDEF > 

IF  (Q)  GOTO  100 
STOP  i 123456 
C 

980  URI TE( 1, 1030) 

1030  FORMAT( '  ILLEGAL  MK '  ) 


o  o  o 


C  PROGRAM  MAP  FAB 


COTO  5 
END 

C  PROGRAM  MAP  HDRED 


C  PROGRAM  MAP  HOREO 

C  TO  EDIT  THE  HEADER  RECORD  OP  A  MAP  FILE 
LOGICAL  0/ CHANGE 

DOUBLE  PRECISION  D I T EMS< 36 >, DWORD 
DIMENSION  HDRC 70  ),  IPTR< 33 ) 

COMMON /F 1LENC/NPR0GC  4  ) ,  NVERS 

COMMON/HEADER/IHDRC40  ) ,  I P  R  J  C  1  ,  IPRJC2<SCALE<XSCALE<  YSCALE 

*  SCLLNG. SCLLAT, RADEA, BADEB, IHDUHBc  24  ), 

*  NL IN , NPT, NSTEPS , NSTEP2, MNFC , HD  1FFC ,  IHDUMC<  14  > , 

*  CNRC  2,4  >, CCNRC  2,4  >, IHDUMDC6  8),  IHOUMEC  82  4  ) 

COMMON/ INC  ON/ INDUMC 1 0  27  ) 

IINSERT  SYSCOM >A*KEYS 

EQUIVALENCE  C  IHDRC  1 > , HDRC  1  )  ) 

DATA  DITEMS  /8HIPRJC1  -8HIPRJC2  .BHNLIN  -8HNPT 

*  8HNSTEPS  -8HNSTEPZ  , 8HMNFC  -8HNDIFFC  ,8HSCALE 

*  8HXSC ALE  ,8 HY SCALE  , 8HSCLLHG  -8HSCLLAT  >  8  HR  AD  E  A 

*  8HRADEB  ,  8  HC  N R«.  1 ,  1  ) ,  8HC NR <  2  ,  1  >  ,  8  H  CN  R<  1  ,  2  ) ,  8  HC  NR  (  2 ,  2  )  , 

*  8HCN  R  c  1  ,3  )  , 8HC  NRc  2,  3  )■  8HCNR<  1  ,  4  ) ,  8HCHRC  2  -  4  >, 8HGCNC  1,  1), 

*  8H GC  N  C  2  , 1  )  ,  8  HG  C  NC  1 , 2  ) , 8HG CN C 2 , 2  > , 8 H GC NC  1-3>,8HGCNC2,3), 

*  BHGCNU -4  >,8HCCN<  2- 4  >,  8HIPR0J  -8HIDESCR  ,  8  HH  EL  P 

*  8HEND  -8H  / 

DATA  IPTR/ 41-42-81, 82, 83-84, 85, 86, 22, 23-24, 25, 26 -27, 28, 

*  51,  32, 33, 54, 53, 56, 57, 58, 59, 60, El, 62, 63-64, 65, 66, 1,7/ 
DATA  NPR0G/8HMP  HDRED/, NVERS/l/ 


URI TEC  1  -  1000) 

1B00  FORMAT <  ' MAP  HDRED  9-JAN-S0') 

IBB  CHANGE3 . FALSE 

CALL  N0PEN8  /*OPEN  INPUT  FILE  FOR  REAO/URI TE 
CALL  NHDIN  /*RE AO  /HEADER/  FROM  FIRST  FILE  BLOCK 
Q»YSNO*AC ' HEADER' , 6, A4NDEF  ) 

IF  <  NOT  13)  GOTO  200 

URI TEC  1  ,  1002  X IHDRC J  ), J  =  t , 22  )  /*IPROJ  AND  IDESCP 
1002  FORMATC ' PROJECT  :  '-6A2,';  DESCRIPTION:  '-16A21 

J  A" -3  9 
8B  JA-JA+40 

IF  < IHOUMEC JA  )  EQ  0  )  GOTO  200 
J8-JA+28 

UR  I  TEC  1 , 1003 >C IHDUMEC J >, J  =  JA , JB >  / • UP  DA  T  E  LINE 

IBB 3  F0RMAT<4A2,I2,l3,lX.A2,'i',A2,'i'A2.I5,2X,19A2> 

GOTO  80 

200  URI TEC  1  ,  1008  > 

108  8  FORMATc  '  ITEM:  '  ) 

READc  1,1102  1DU0RD 
1 102  FORMA  TC  A8  ) 

DO  210  J»l , 36 
I TNO-J 

IF  CDUORD  EQ  DITEMSC  J  ))  GOTO  220 
210  CONTINUE 

URITEC I, 1002) 

10B7  FORMATC' ITEM  NOT  RECOGNISED  TYPE  “HELP"') 

GOTO  200 

220  IF  CITNO  EQ  33  OR  ITNO  EQ  36)  GOTO  600  / • E N D  OR  <CR> 

IF  CITNO  EQ  34)  GOTO  500  /•"HELP" 

!SU»l  /*D I TEHS  1  TO  8  -  INTEGERS 

IF  CITNO  GT  8)  ISU»2  /*DIT£MS  9  TO  31  -  REALS 

IF  CITNO  GT  31)  ISU-3  /»DITEMb  32,  33  CHARACTERS 


.  ■  S.  N. 


O  Cl 


C  PROGRAM  MAP . HDRED 


MS" IPTR< ITNO) 

COTO  <  221. 222,223), I SU 

221  URITEC  1,  1221  )OITEMS<  ITNO).  IHDRf  NS  ) 

1221  FORMATC A8. '=',118) 

COTO  238 

222  URITE<  1,  1222)OITEMS<  I THO ) , HDRC HS  ) 

1222  F0RNATCA8, '  = ' ,  F  28  10) 

GOTO  238 

223  IF  (ITNO. EQ. 32)  NE  =  6  /"END  WORD  OF  IPROJ 
IF  (ITNO  EQ  33)  NE=22  /"IDESCR 

UR  I  TEC  1,  1223)01  TEMS(  1 TNO > , C I HDRC J  ), J  =NS, ME  ) 

1223  FORMATC A8-  ' «' , 16A2  ) 

230  0*YSN0«AC ' CHANCE  I T '  , 9 , A* N OE F  ) 

IF  (  NOT  .Q  >  COTO  208 
CHANCE"  TRUE 

URI TEC  1,  1088) 

1008  FORMAT < ' NEW  VALUE:') 

COTO  <  231, 232,233), ISU 

231  READ<  1,1231  HHDRCNS) 

1231  FORNATC110) 

GOTO  208 

232  READ(  1 ,1232  )HDR ( NS  ) 

1232  FORMA  TC  F20 . 2 ) 

GOTO  288 

23  3  REAO<  1  ,  123  3  )(  IHDRC J  > , J*NS , NE > 

1233  FORMA T ( 1 6A2  ) 

GOTO  280 

C 

300  URITEC 1, 1508 ) D I T EMS 
1300  FORMATC  '  ITEMS  AVAILABLE  ARE  s  -  '  ,'8(  2X  ,  A8  )  ) 

GOTO  200 
C 

600  1F(  NOT. CHANCE)  GOTO  610 

CALL  UPDATE  /"ONLY  UPDATE  IF  A  CHANGE  HAS  EE  E  N  MADE 
CALL  NHD8AK  /*URITE  HEADER  BACK  TO  FILE,  IF  CHANGED 
610  CALL  NCLOSI 

O-YSNOIAT ' MORE  FILES’, 10, AINOEF) 

IF  (0)  COTO  100 
STOP  1123436 
END 

SUBROUTINE  NOPENB  /"OPEN  A  MAP  FILE  FOR  BINARY  INPUT  OR  OUTPUT 
DIMENSION  IFILEIC2Q) 

LOCICAL  Q 

COMNON/INCOH/INOEV, INPTR.  INSZ,  INBUFi  1024  ) 

•  INSERT  S  YSC  OM  > A  IKE  YS 


I NDEV  =  5 
INSZ" 1024 
NLEN-40 

108  0"0PNP*A< '  INPUT  FILE'  ,1B.  AIRDUR, IFILEI ,NLEN, INDEV) 

IF  (  . HOT  Q  )  GOTO  108 

RETURN 

END 

C* •**•••••»••**••#»*»»*»••• •**»•*»*•**»••*«**.***»•*** ••••••••« 

SUBROUTINE  NHDBAK  /«URITES  THE  FIRST  BLOCK  INTO  /HEADER/ 
COMMON/HEADER/I HDRC  1024  > 

COMMON /I NCON/ INDEV- INPTR.  INSZ,  I NBUFc  1024  ) 

•  IMSERT  SYSCOM  >KE YS . F 


ono 
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C  PROGRAM  HAP.HORED  i 

I 

i 


f 


DO  1B0  J*1 . INS 2 
IBB  INBUP( J  )-IHDR< J  ) 

CALL  PRUF* t< K*UR  IT+K$PREA ,  INDEV,L0C<  INBUF  ),INS2.INTL<B). 
•  IRNU i ICODE  ) 

IF  (ICODE. EQB  AND  1RNW.EQ  INS2)  RETURN 
URITE< I. 1B2BHRNU, ICODE 
1  B2B  FORHAT<  '  PRWFM  FAILURE  IN  NHDBAK  :  '  .  2  18  > 


n  o  p» 
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C  PROGRAM  map  hplot 


C  PROGRAM  MAP  HPLOT 

C  TO  PLOT  A  MAP  FILE  OH  A  MEULE T T - F A C K A RD  PLOTTER 
LOGICAL  Q-QC 
INTEGER**  DPOS 

CONMOH/COORDS/’MK  -X  .  V,  NPTS,  NPT2.  NFC.  FCODE(  1  7  > 
COMMON/HEADER/IMDUMl  4B)-IPRJC1-1HDUMJ, 

*  SCALE . IHOUMA( 36  >-CNR(2-4  )  -  I  HD  U  MB  (  9  B  8  ) 
IINSERT  SYSCOM>A*KEYS 


UR  I TE  <  1.  1  BBB  > 

1  BOB  FORMAT< ' MAP  HPLOT  22-JAN-8B') 

IBB  CALL  HPIHIT  /*INITIAL12E  PLOTTER 
1 1 B  UR  I TE  < 1.  1 BB  2  > 

1BB2  F ORMAT ( ' PEN  NUMBER:') 

READ<  l , 1  IBB  >J 
1  IBB  F  ORMAT  t 16) 

IF  (J  LT.l  OR  J  CT.4)  GOTO  11B 
CALL  HPEN(J)  /'•SELECT  PEN  J 
116  Q*RNUM«A< ' CROSS  SIZE  <  I  6  )  '  .  1 5, A*DEC . DPOS  ) 

IF  ( . NOT  0  )  GOTO  1 16 
NCROSS«INTS(DPOS ) 

CALL  HOPENI  /*OP£N  INPUT  FILE 
CALL  NHOIN  /  *  I  N  PU  T  THE  HEADER  BLOCK 

131  U»CHR( 1, 3)-CNR(  1  - 1 >  Z*UIDTH  OF  FRAME 
H  =  CNR( 2. 3)-CNR(  2 , 1 )  /'•HEIGHT  OF  FRAME 
IF  (H.LE  B  OR  U  LIB)  COTO  950 
QC*YSNO«A( 'SPECIFY  CORNERS'- 15-A*DN0  ) 

IF  (  HOT  .80  GOTO  6B 

132  Q*RHUHtA('S  U  X ' , 6 . A tOEC , DP  OS > 

IF  < . NOT . Q  )  GOTO  152 

IF  (DPOS  LT.B  OR  DPOS  CT. 3080  )  GOTO  152 
SUX«FLOAT<  DPOS) 

133  0»RNUM*A('S.U.  Y ' , 6, A4DEC , DPOS  ) 

IF  < -  NOT . Q  )  COTO  153 

IF  (DPOS  LT  .B  OR  DPOS  CT  2080  )  GOTO  153 
SUY-FLOAT( OPOS) 

154  Q»RNUM*A(‘N  E  X ' . 6 . AIDEC . DP  OS  ) 

IF  ( . NOT . Q  )  GOTO  154 

IF  (DPOS  LE  IFIX(SUX)  OR.  DPOS  GT  30BB >  GOTO  154 
RNEX«FLOAT(  OPOS  > 

153  Q*RHUM«A('N  £  Y  '  , 6 , A «DEC . DPOS  ) 

IF  ( . HOT  0 )  COTO  155 

IF  (DPOS  LE  IFIX(SUY)  OR.  DPOS  GT  2B0O)  GOTO  155 
RNEY=FLOAT (  DPOS ) 

PU-RNEX-SUX 
PH*RNEY-SUY 
COTO  70 

SB  IF  (SCALE  LT  IB  0)  COTO  94B 

39  Q»RNUM*A< ' OUTPUT  SCALE  ( I  9  )'  .  1 7  .  A*DEC . DPOS  ) 

IF  ( .NOT  0 )  GOTO  59 

IF  (DPOS  LT  IBB  OR  DPOS  CT  lOBOBOBBO)  GOTO  39 
XSC«SCALE/'FLOAT<  DPOS  )«B  .  D  -'*PLOTTER  UNITS  ARE  1 /'8  MM 
IF  ( U * XSC . G T  3BBB  B  OR  H*XSC  CT  2BB0  0>  GOTO  93B 
GOTO  79 

6B  SUX*0  B 
SUY-B  B 
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C  PROGRAM  MAP  HPLOT 


IF  <IPRJC1.ES. 99)  COTO  65  /-IMACE-FILE  UMITS 
0-YSN04AC ' SPECIFY  SCALE' . 13. A*MDEF> 

IF  <S)  Coro  58 
65  PU-3B00  0 
PH-2000 . 0 

70  XSC-PU/U  /-X  SCAL1HC  FACTOR  TO  FIT  PLOTTER 
YSC-PH/H  /  *Y  SCALIHG  FACTOR 
IF  (OC)  COTO  80 

XSC=AMIN0< XSC.YSC)  /♦CONVERSION  FROM  FILE  (MM)  TO  TEKTRONIX  UNITS 
IF  (SCALE.LT  10.0)  COTO  79 
OSCALE-SCALE-8. 0/XSC 
URI TEC  1.  1004  )OSCALE 
1004  FORMATC ' SCALE  PLOT TE D= ' , F 1 6  2 > 

79  YSC-XSC 

80  0-YSN0«A( ■ FRAME'  .5 . A4NDEF  ) 

IF  < . NOT  B  )  COTO  95 

CALL  NPMOVE 

DO  90  J»l,5  /-DRAU  THE  FRAME 
JA«  J 

IF  < JA  EO. 5  )  JA-1 

IX* IFIXC  CNRC 1 . JA  )-XSC-0  5-SWX) 

IY-IFIXC  CNRC2.JA  )-YSC-B . 5-SUY  ) 

CALL  HPOSPT < IX.  I Y  > 

90  CONTINUE 

95  ITEM-0  /-COUNTS  ITEMS  (LINE,  POINT)  IN  A  SHEET 
2B0  CALL  NFETCH 

IF  CMK  CE  6  OR  HK  LT  0)  GOTO  600 

IF  ( MK  Eg  3  OR  MK.ea  4  OR.  MK.EO  5)  ITEH-ITEM-1 
IXP-IFIXCX-XSC-0  5-SUX) 

IYP-IF1X(Y*YSC»0 5-SUY  ) 

COTO  (250, 250, 300, 400, 400), MK  /-DUMMY  FOR  MK-2 
C  HK -0  OR  NK«t  INTERIOR  POINT  OF  A  LINE  OR  END  OF  LINE 

250  IF  (IXP-IXPO.EO  8  AND  IYP-IYPO  EQ  0)  COTO  2BB  /-SAME  POINT 
CALL  HPOSP  T ( IXP,  1YP)  /-DRAU  TO  X.Y 
GOTO  310 

C  HK -3  START  OF  LINE 

380  CALL  HPMOVE  /-SET  VECTOR 

CALL  HPOSPT( IXP,  IYP)  /-MOVE  TO  X.Y 
310  IXPO-IXP  /-PREVIOUS  POINT 
IYPO-IYP 
COTO  200 

C  MK -4  OR  5  POINT  OR  FIDUCIAL 

400  CALL  PLOTCRC IXP, IYP. NCROSS >  /-PLOT  POINT  AS  A  CROSS  (MY  SUBR  ) 

GOTO  208 

C  MK -6 .  END  OF  THE  MAP 
600  URITEC  1.  1 00  3  )MK ,  ITEM 
1003  FORHATC ' MK-' ,  16,  '  AFTER  ITEM',  16) 

6SB  CALL  HPEND  /-PURGE  BUFFER,  IE  DRAW  EVERYTHING 
CALL  NCLOSI  /-CLOSE  INPUT  FILE 
0»YSN0*A< ' ANY  MORE  F I LE S ' . 1 4 , A t NDEF  ) 

IF  <Q>  COTO  188 
STOP  i 123456 
C 

930  URI TEC  1, 1022) 

1022  FORHATC'HAP  TOO  LARGE  AT  THIS  SCALE') 

GOTO  59 
C 

940  URITEC 1, 1020) 

1020  FORHATC ' INPUT  SCALE  MOT  SPECIFIED') 

GOTO  151 


I 


^  N. 


44 


C  PROGRAM  MAP . HPLOT 


C 

930  UR  I  TE(  1 ,  1021  >CNR<  1.1  )  .  CNR<  2.  1  ) ,  CNR(  1  .3  >,  CNR*.  2  .  3  > 

1021  FORMAT ( '  FIDUCIALS  I MCORRECT ' , 2* /2F 1 0  2  >  ) 

GOTO  630 
END 

t************** . .  . . . 

SUBROUTINE  PLOTCR< IX, IY,N >  /*PLOTS  A  POINT  AS  A  CROSS,  FOR  FIPLOT 
C 
C 

CALL  HPMOVE 

CALL  HPOSPTC I X+N , I Y ) 

CALL  HPOSPTt IX-N, I Y> 

CALL  HPMOVE 

CALL  HPOSPTC IX, I  Y»N  > 

CALL  HPOSPTt IX. I Y-N) 

RETURN 

ENO 


►  * 
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C  PROGRAM  MAP  IN 


C  PROGRAM  MAP  IN 

C  TO  READ  AN  “EXCHANGE  FORMAT"  MAGTAPE  INTO  A  MAP  FILE 
LOGICAL  NOCORN 
DIMENSION  !BUF2(2O0> 

COMMON/COORDS/MK,X,Y, NPTS.NPT2, NFC. FCODE!  1?  ) 
COMMON/FILENC/NPROG!  4  >.  NVERS 

COMMOM/HEADER/IHDR! 1O0),CNR( 2, 4  >,  IHDUMh!  908  > 
COMMON/INTRDC/IBP,IBUFC20B> 

CONNON/OUTCOM/IOTDUM(  182?) 

DATA  NPR0G/8HNAP  IN  /.NVERS/I/ 


ONIT-B  02  /•CONVERTS  EXCHANGE  UNITS  TO  MM 
WRITE! 1, 1000) 

1000  FORMAT!'  MAP.  IN  9-JAN-30',"  TAPE  UNIT  NO:') 

READ!  1,1100  )MTUN  1  T 
1  100  FORMAT! 16) 

CALL  REN! MTUN  IT  ) 

NBLOCK»0  /•COUNTS  BLOCKS  ON  MAGTAPE 
NF I LE  =  0  /*COUNTS  FILES  FROM  THIS  MT 
C  START  NEXT  OUTPUT  FILE 

98  LIN£  =  0  / *COUN  TS  BLOCKS  WITHIN  ONE  HaP 

I TE M“ 0  /•COUNTS  ITEMS  CLINE,  POINT)  IN  A  SHEET 
NOCORN3 . FALSE  .  /*FOR  FIDUCIAL  CHECK 
NFILE-NFILE+1 

100  I ND  =  6  0  /*  EXPECTED  HO  OF  WORDS  !120  CHAPS) 

CALL  MTREAOCHTUNIT  ,  NB  LO  CK  ,  IND)  ,'*INCREMENTS  NBLOCK  AS  WELL 
LINE-L1NE+1 

IF  CIHD  E0.4)  STOP  : 1 000  /*EOT 
IF  (IND  Et  3)  GOTO  990  /*EOF 
IF  CLINE  GT  .  10)  GOTO  110 

IF  (LINE  HE  1  OR  NBLOCK  ED  l>  GOTO  102 
MK= INT  RD<  2  ) 

IF  CMK  EO  30)  GOTO  980 
WRITE! 1, 1006) 

10B6  FORMAT!'  ANOTHER  FILE  ON  TAPE  TYPE  S  TO  CONTINUE') 

PAUSE  /*T YPE  *S“  TO  CONTINUE 

102  DO  103  J»l,120  /*PREPARE  IBUF  FOR  PRINTING 

103  1BUF2C J )«LS! IBUF! J  ). 8  )  /‘SHIFT  INTO  CORRECT  BYTE 
WRITE!  1,  1BB2)!IBUF2(J  ),J=1,120)  /*t>0-UORD  BLOCK 

1002  FORMAT!  12BA1 ) 

IF  CLINE. HE. 10)  GOTO  100 

CALL  HOPENO  /* OPEN  THE  OUTPUT  FILE 

IF  CNF1LE.E0  1)  CALL  PERMHD  / *GE  T  PERMANENT  HEADER  DETAILS 
DO  108  J *4 1 ,1024 

108  I HDR!  J  )»B  /♦CLEAR  HEADER  EXCEPT  FOR  PERN  DETAILS 
CALL  UPDATE  /• DETAILS  OF  THIS  RUN 
GOTO  100 
110  T  TEH* I  TEH* 1 
HK-INTR0C2  > 

IF  <HK  CT  6  >  GOTO  9B0 

IF  CITEM  LE.4  AND  MK.NE  3)  GOTO  B?3 
IF  CITEM  CT.4  AND.  MK . EO . 5  )  GOTO  880 
106  J«INTRD! 2)  /«SHOULD  BE  ZERO 
IF  (J.NE  0)  WRITE! 1, 1003)0 

1003  FORMAT!'  IC0R-',16,'  IN  ITEM', 16) 

IX- INIRDC6) 

IY-INTRD! 6 ) 

X- FLOAT! IX ) *UNI T 
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C  PROGRAM  MAP. IN 


Y«FLOAT( IV ) *UNI T 
NFC-INTRDC 2  ) 

IF  (NFC  LT . 1  )  GOTO  130 
IF  (NFC.CT. 17)  NFC-17 
DO  126  J-l.NFC 
JA-INTR D(6  > 

126  FCODE(J)*FLOAT(JA) 

130  IF  (1TEH.LE.4)  URI T£<  1 .  1004 > I T E N , MK ,  1 X , 1 Y  ,  NFC 
1004  F0RMAT(I6, '  >'  .14,218,  14) 

1008  FORNAT< 16. '  i ' , I4,2F8  0. 14  ) 

GOTO  (  900,  900,300, 400,500, 600), MK 

300  CALL  NSTORE  /‘MK-3  ITEM 
NK-0  /‘INTERIOR  LINE  POINT 
NPTS-1  /‘FIRST  POINT  WITH  HK-3 
NPT2-0 
NFC-0 

310  IND-60 

CALL  MTR£AD(HTUNIT , NBLOCK ,  I N  D  > 

IF  (IND.NE  1)  NRITEC 1,1005>IND, NBLOCK 
1003  F  ORMAT  < '  I  NO- ' , I  6 , '  IN  BLOCK*.  16) 

DO  320  3-1,30 
I  OX- 1 NTRO< 2  >-30 

IF  (IDX  EO.-50)  GOTO  340  /‘END  OF  LINE 
I DY-I NTRD< 2  >-30 
IF  (IOY.EB.-30)  GOTO  340 
NPTS-NPTS‘1 

IF  (NPTS.LT. 10000)  GOTO  316 

NPTS-0 

NPT2-NPT2M 

316  IF  ( NPTS . EQ  .  2  AND  NPT2  EQ  0)  COTO  318 
X-XO 
Y-YO 

CALL  NSTORE  /‘PREVIOUS  POINT,  MK-0 
318  XO-X‘FLOAT(  IDXXUNIT 
YO- Y‘FLQAT (  IDYXUNIT 
320  CONTINUE 
COTO  310 

340  IF  (NPTS  EQ  1  AND  NPT2  EQ  0)  GOTO  35B 
MK- 1  /‘END  OF  LINE 
X-XO 
Y-YO 

CALL  NSTORE  /‘END  POINT  OF  LINE 
COTO  100 

330  URI TE(  1,  1BB7)ITEH 

1007  FORMAT( ' ITEM' , 16, '  IS  A  ONE-POINT  LINE') 

GOTO  100 

400  MK- 4 

CALL  NSTORE 
GOTO  IBS 

580  CNR(  1 , ITEM  ) -X 

CNR( 2 , 1  TEN )-Y  /‘FIDUCIALS 
IF  (NFC  ME  0)  «RITE(  l , 1089  )I TEM-NFC 
1 0B9  FORMAT!'  ITEM'  ,16,'  HAS  NFC-',  14) 

IF  (  ITEM  LT  .4  )  GOTO  IBB 
IF  (NOCORN)  GOTO  930 

CALL  HSTART  /‘PLACE  /HEADER/  INCL  CNRS 
HK-3  /‘FIDUCIAL 
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C  PROGRAM  hap.  in 


NFC-8 

DO  528  J-i >4 
X  =  CNR<  1 ,  J  > 

Y»CNR( 2, J  ) 

520  CALL  NSTORE 

coto  tea 

c 

680  HK-6 

CALL  NSTORE 
CALL  NCLOSO 
URI TE(  1.  1020  )NK .  ITEM 
1  020  F ORMAT (' MK . 16/  '  AFTER  ITEM'.  16) 

GOTO  98 
C 

878  NOCORN-  TRUE 
GOTO  80S 

sea  MK-4  7*MK  5  IN  BODY  OF  DATA,  CHANCED  TO  MK  4 
865  URI TE( 1, 1B98>NK. ITEM 
1098  FORMATC'  WARNING  MK-',13.'  IN  ITER',16) 

GOTO  106 

900  WRITE( 1, 1098>NK, ITEM 
GOTO  100 
C 

950  URI TE<  1,  1021  ) 

1021  FORMATS  F I DUC I ALS  INCOMPLETE') 

CALL  NCLOSO  /-USE  NDELET  WHEN  AVAILABLE 
COTO  990 

980  WRITE( 1, 1822>MK, NBLOCK 

1022  FORMATS  MK--.I3,'  IN  BLOCK',  16) 

990  CALL  REUCHTUNIT  ) 

STOP  i 123436 
END 

£•*«•••*»**»*»*••*»••**»*»•*•*••***»•*******»***«*»*»•**•**••»*»••***»**•••»** 

SUBROUTINE  HTRE AD( NT  UNI T . NBLOCK . I ND  ) 

C  SUBR  TO  READ  ONE  BLOCK  FROM  MT  AND  PLACE  EACH  CHAR  INTO  A  UORD  OF  IBUF 
C  IND-l:  OK  2  >  NU  WRONG  3>  EOF  4>  EOT  5<  ERR 
LOGICAL  EOF, EOT, ERR 
COMMON / 1  NT  ROC/ 1  BP,  IBUF(  1  > 

C 

C 

CALL  MTRE<MTUNIT, IBUF, 200, EOT, EOF, ERR. NW>  /•HTIOLJB  SUBR 
NBLOCK-MBLOCK+1 
IF  (NOT  EOF)  GOTO  18 
WRI TEC  1,  1000  )NBLOCK 
1 0B0  F  ORMAT ( '  EOF  IN  BLOCK', 16) 

IND-3 

RETURN 

18  IF  i  NOT  EOT)  GOTO  28 
WRI TE( 1, 1BB1  )NBLOCK 
1BB1  FORMAT ( '  EOT  IN  BLOCK',  16) 

I  NO  -4 
RETURN 

28  IF  (NOT  ERR)  COTO  30 
WRI TE(  1,  1002  )N6 L  OC K 

1882  FORMAT ( ’  ERROR  IN  BLOCK'. 16) 

IND-5 

RETURN 

3B  IF  (NW  EQ  IND  >  GOTO  48 
WRI TEC  1, 1BB3)NW, NBLOCK 

1883  F ORMA T ( 16, '  WORDS  IN  BLOCK', 16) 


C  PROGRAM  MAP.  IN 


imo-2 

RETURN 

C 

48  00  38  J-1,NW  /‘UNPACK  CHARS  TO  WORDS »  BACKWARDS 

K-NU-J‘l 
KA-K‘K 
N-I8UFCK > 

N1-RT(N.?>  /‘RIGHT  TRUNCATE  OR  USE  ;  AMD<N,12?> 

H  2  *  R  S  ( N, 8)  /‘RIGHT  SHIFT  8  PLACES 
N2-RT (N2,7> 

I BUF< KA-I  )«N2 
1 8UF< KA)-N1 
38  CONTINUE 
IND-1 

I  BP- 1  /‘READY  FOR  USE 

RETURN 

END 

C* 

FUNCTION  INTRD(N) 

C  CONVERTS  NEXT  N  <1  TO  6)  CHARS  TO  AN  INTEGER 
LOGICAL  NUHST -NEC. ERROR 
DIMENSION  I STRNGl 6  ) 

COMMON/INTRDC/IBP,  IBUF(  1  > 

C 

c 

ERROR- .FALSE. 

IF  (N.LT.l  OR.  NOT  8)  RETURN 

HUH-B 

HUHST-  FALSE.  /‘NUMBER  NOT  STARTED 
NEC-FALSE.  /‘HO  NIHUS  SIGN  SEEN  YET 
DO  188  J-l.N 

NEXCH«RT( IBUF( IBP‘J-1 >, 7>  /‘LAST  7  BITS  ONLY 

IF  ( NEXCH  EQ  .  32  OR  NEXCH  EQ  43  OR  NEXCH  EQ  45)  GOTO  SB 

NEXCH-NEXCH-4B  /‘CONVERT  TO  NUMBER 

IF  CNEXCH.LT. 8  OR  NEXCH. GT. 9)  COTO  900 

IF  (NUN  CT. 3276)  GOTO  9BB  /‘OUTPUT  MUST  NOT  EXCEED  3276B 
HUM-HUM* 18‘ NEXCH 
COTO  68 

38  IF  ( HUHST )  GOTO  9BB 

IF  (NEXCH. EQ  43)  GOTO  68  /*  PLUS  SIGH 

IF  (NEXCH. HE. 43 >  GOTO  188  /*  MINUS  SICH 
MEG-TRUE 
68  HUHST-  TRUE  . 

188  CONTINUE 

IF  (HEC)  NUN--NUN 
I NTRD-NUM 
138  I  BP- 1 BP+N 
RETURN 
C 

9BB  DO  918  J-l ,N 

918  ISTRHGC J )-LS( IBUF( I8P‘J-i >,8 >  /‘MOVE  TO  LEFT  BYTE 
URI TEC  1 < 18B8)( ISTRNG( J >. J-l.N) 

1  BBS  FORMAT( ' ILLEGAL  INTRD  STRING;  '.6A1> 

IHTRD-8 
ERROR-  TRUE  . 

GOTO  138 
END 


/*SP  ♦  - 
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C  PROGRAM  MAP. LIST 


C  PROGRAM  MAP. LIST 
C  TO  LIST  A  MAP  PILE 

LOGICAL  0- IMTPTS 

COMMON/HEAOER/IHORC 100>,CNR( 2,  4  >,  GCNRC  2,  4  >, 1 HDUMAC  68 > , I HDUNBC  824  > 
C0HM0N/INC0M/1NDUHC 1B2? >  /*NOT  NECESSARY? 

COMMON/COORDS/HK ,X, Y, NPTS. NPTZ, NFC. FCODEC  1?  > 

*  INSERT  SYSCOM  >A$KE YS 
C 
C 
C 

URITEC 1. 1800) 

1080  FORMAT< ' MAP  LIST  9-JAN-80'  > 

IB  CALL  NOPENI  /*OPEN  INPUT  FILE 

CALL  NHOIN  /*R£AO  /HEADER/  FROM  FIRST  FILE  BLOCK 
0=>YSN0*AC  '  HEADER  6,  ARNDEF  > 

IF  < . NOT . Q  )  GOTO  9B 

WRI TE< 1 < IB  1 2  >( I HDRC J  )- J-l . 22  > 

1012  FORMATC ' PROJECT :  '-6A2-'  DESCRIPTION!  '-16A2) 

JA--39 

8B  J  A“ JA ♦ 4B 

IF  < IHDUMBC JA  )  EQ. 0  )  GOTO  82 
J  B“ J A*  28 

URITE<  1,  101 1 )< IH0UM8C JJ-J-JA-JB)  /*  UPDA  TE  LINE 
1011  FORMATC 4 A2, I2-I3-1X-A2- ' > ' -A2- '  s'A2- I5-2X-19A2  ) 

COTO  80 

82  WRI  TEC  1.  1010MCNRC  1  ,  J  ),  CNRC2- J  >,  J-l  ,  4  ) 

1010  FORMATC ' CN R S '  -8F9  2) 

WRI TEC  1-  1013)CGCNRC1-J),GCNRC2- J),J»1,4> 

1013  FORMATC ' GCNRS'. 8F9 .3  ) 

90  Q*YSH04AC'DATA' -4-AfNOEF) 

IF  (.HOT  II)  GOTO  900 

I MTPT S»YSNO*A< ' INTERIOR  LINE  PO  I N TS ' , 2 0 . A «NDEF > 

MK-3 

NLIN«8  /*N0  OF  LINES  IN  THE  FILE 
NPT  =*0  /  »N0  OF  POINTS  IN  THE  FILE 

100  MKO*MK 

CALL  NFETCH 

GOTO  C 110- 120-125- 130-140. 150- 160), MK+1 
110  IF  CMKO  EO  0  .OR  MKO.EO  3)  COTO  112 
URITEC  1. 1020  )MK . MKO 
1020  FORMATC ' HK 12,  '  AFTER  MK '  ,  I  2  ) 

PAUSE  Ull 

112  IF  CINTPTS)  URITEC  1.  1003)MK,X,Y 
116  NPT  SA  =  NP  TS A  +  1 

IF  CNPTSA  LE  9999)  GOTO  100 

NPTZA»NPTZA+1 

NPTSA=B 

GOTO  100 

120  IF  illk'O  EO  0  OR  MKO  EO  3>  GOTO  122 
URITEC  1 .  1020  >HK . MKO 
PhUSE  mi 

122  URITEi 1 . 1003 )MK , X  Y, NPTZA . NPTSA 
GOTO  116 

123  WRI TE<  1  -  1020  )MK - MKO 
PhUSE  nil 

GOTO  100 

130  IF  ,  MKO  EO  1  OR  MKO  EQ  4  OR  MKO  EO  5>  GOTO  132 
URITEC  1 ,  1020  )MK . MKO 
PAUSE  nn 
132  NLIN-NLINM 
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C  PROGRAM  MAP. LIST 


IF  (NFC.LE  0)  URITE(  1 , l 003  INK . X , Y ,  HPT?,  MPTS,  MFC 

IF  (MFC  GT  8  >  URITE(  1  ,  1003  >MK,  X  .Y  ,MPTZ  UPTS,HFC,  '  FCOl'E(  J  ><  J=  1  NFC 

1003  FORMAT(I2,2F10.2,I3,  14,  13.  17F1Q  1) 

MPT  SA  =  2  / *  TH  I S  POIMT  AMO  END  POINT  (NK=1> 

MPT  ZA  =  0 
GOTO  100 

MB  IF  (MKO  EO.I  OR.  MKO.EO  4  OR.  MKO  EO  5>  GOTO  142 
URITE(  1,  1020  >MK , MKO 
PAUSE  : 1 1 1 

M2  IF  (NFC  LE  0)  WR  1  TE(  1  ,  1  004  )HK ,  X  ,  V  -  MF  C 

IF  (MFC  CT.0)  MRITE<  1  -  1004  >MK, X,Y, MFC. < FCODE( J  >, J= 1 -MFC  ) 

1004  F  ORMAT  (  12, 2F10  2,7X. 13, 17F10  1 > 

GOTO  154 

150  IF  (MKO  E0.5)  COTO  152 
URI TE(  1 ,  1020  >MK , MKO 
PAUSE  <111 

152  URI T E (  1,  100  3 )MK ,  X  •  Y 
154  NPT-NPTM 
GOTO  100 

160  IF  (MKO  EO.I  OR  MKO  EQ  4  OR  MKO  EQ  5)  COTO  162 
URI TEC  1 ,  1020  )MK , MKO 
PAUSE  :  1 1 1 
162  URI TE(  1 ,  1  003  )MK 

UR  I T  E  (  1,  1805TNPT ,MLIN 

1005  F  ORMAT  ('FILE  CO N T A  I H S ' ,  I  6 ,  '  PO I N T S .  '  , 1 b  ,  '  LINES  ') 

900  CALL  N  CL  OS  I 

Q=YSNO*A( ' MORE  F  1  L ES '  ,  1 0 ,  A *N DE F  ) 

IF  (Q  >  GOTO  10 
STOP  :  123456 
END 


< 

i 
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C  PROGRAM  NAP  OUT 


C  PROGRAM  NAP  OUT 

C  TO  OUTPUT  A  NAP  FILE  TO  A  MAGTAPE  IN  "EXCHANGE  FORMAT" 
LOGICAL  ERROR , Q 

COMMON,' 'COORDS/MK.X,  V,  NPTS.  HP  T2 ,  NFC,  FCODEf  17  ) 
COMMON, 'ERR/ERRQR 
COMMON/HEADER/IHDR!  1 0  24  > 

COMMON, 'INCOM/INDUM! 1 02?  ) 

C0MM0N/INTRDC/I8P,  I8UF!  20  0  ) 

IINSERT  S YSCOM  >A*KE YS 
C 
C 
C 

U NI T  =  50  0  /*COHVERT  MM  TO  EXCHANGE  UNITS 
UNITE! 1, 1800) 

1000  FORMAT!'  MAP  OUT  9-JAN-80'."  TAPE  UNIT  HO:') 

READ!  1,1100  )MTUN  1  T 

1  100  FORMAT! I  6  ) 

CALL  REUCMTUNIT)  /*MTIOLIB  CALL 
NBLOCK=0  * C 0 U N T S  BLOCKS  ON  MAGTAPE 
100  CALL  NOPENI  ,'*OP£N  INTERNAL  FILE 
CALL  NHDIN  /'♦READS  FIRST  BLOCK 
C 

WRITE!  1,1001) 

1001  FORMAT! ' UR  I TE  TEN  HEADER  LINES:') 

DO  110  J  =  1 -  10 

READ!  1,1101  )!  IBUF!  1),  1  =  1,  120) 

1101  FORMAT! 1 20A  1  ) 

DO  104  K  =  1 ,  120 

104  IBUF! K  )  =  RS!  IBUF! K),8  )  /*SHIFT  TO  RICHT 

CALL  MTURI T! MTUN IT , NBLOCK >  /*MTIOLIB  SUBROUTINE 
IF  (ERROR)  GOTO  993 
110  CONTINUE 
C 

200  DO  201  J  =  1 ,  120 

201  IBUF<  J  )=  :240  *  F  I  L  L  IBUF  WITH  SPACES 

I  BP  =  1 

202  CALL  NFETCH  /*GET  NEXT  ITEM 
MKIND=MK*1 

GOTO  !  210,  210,210, 240,240, 240, 260  ), MKIND 
C 

210  X  =X  «U N  1  T  +  0  5 

IF  !  X  LE  32767  0  AND.  X  GE  -32768  )  GOTO  212 
URITE!  1,  10  03  )X 
X=0  0 

212  I X= IF  I  X! X  ) 

Y • Y  *U N  I  T  +0  5 

IF  ! Y  LE  32767.0  AND  Y  GE  -32763)  GOTO  214 
WRITE!  1  ,  1003  )Y 
Y  =  0  .0 

214  I  Y*  IF  I  X!  Y  ) 

IDX-IX-IXO  /'•CHANGE  SINCE  PREVIOUS  POINT 
IDY-IY-IYO 

IXO=IX  /‘SAVE  PRESENT  POINT,  FOR  PREVIOUS 
I  Y  0  =  I  Y 

I AX  =  I ABS! I DX  > 

IAY=IABS! IDY) 

IF  ! I A  X  LT  50  AND  1  AY  LT  50)  GOTO  220 
C  FOR  I DX  OR  IDY  MORE  THAN  49 
N  UM  *M  A  X8 !  I  AX,  I  AY  ) 

NUM=< NUM-I  )/49+ 1  /‘DIVIDE  INTO  NUN  PIECES 
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C  PROGRAM  MAP  OUT 


RNUM=FLO AT ( NUM  > 

DX=FLOAT( I  OX  )/RNUM 
DY=FLOAT( IDY  )/RNUM 
DO  216  J=1 , HUM 

X  “F  LO  A  T< J  )*DX  / «£  XA  C  T  LOCATION 

IX-IFIXC X  +  0  5  )  /♦NEAREST  INTEGER 

IDX=lX-IXO  /*CHANGE  SINCE  PREVIOUS  POINT 

IXO  =  IX  / *N £ W  PREVIOUS  POINT 

Y  aFLO  A  T(  J  )  *  D  Y 

IY»IFIX(  Y*0  5  ) 

IDY-I Y-I YO 
I  YO  =  I  Y 
COTO  224 

216  CONTINUE 

217  IF  <  H  K  £Q  0  )  GOTO  202 

218  CALL  INTURT(B>4  )  /*PACK  REST  OF  BLOCK  WITH  ZEROS 
I  PC  TR= IPCTR  +  1 

IF  (  1PCTR  LT  30  )  COTO  218 
CALL  MTURI T( MTUN I T - NBLOCK  ) 

IF  (ERROR)  GOTO  998 
GOTO  208 
C 

228  IF  (IDX  EQ  0  AND  IDY  EM  AND  MK  EG  0)  GOTO  202 
NUM  =  1 

C  OUTPUT  AN  INTERIOR  LINE  POINT,  MK>=0  OR  1 
224  IDX»IDX*3B 

CALL  INTURTC IDX >2) 

IDY=»IDY  +  30 

CALL  INTURT< IOY, 2) 

IPCTR-IPCTR+1 

IF  ( IPCTR  LT. 30  )  GOTO  228 

CALL  MTURI  HMTUNIT  ,  NBLOCK  ) 

IF  (ERROR)  GOTO  999 
I  BP  =•  1 
I PCTR“0 

220  IF  (NUM  GT  1  )  COTO  216 
COTO  217 
C 

240  CALL  IHTURT(Mk, 2  )  /•POINT,  FIDUCIAL  OR  STmRT  OP  LINE 
CALL  INTURT<0,2)  / • “  1 C  0  R " 

X*X*UNIT‘,0  5 

IF  (X  LE  32767  0  AND  X  GE  -32768)  GOTO  242 

URI TE(  1  ,  1003  )X 

1003  FORMAT( ' COORDINATE  TOO  LARGE  TO  OUTPUT, FID  0  '  UNITS’ 

X*0 

242  IX-IFIX(X) 

CALL  INTURT(IX,6> 

Y«Y*UN  I  T  *0  5 

IF  (Y  LE  32767  0  AND  V  GE  -32768)  GOTO  244 

URI TEC  1  ,  1  0  0  3  )  Y 

Y»0 

244  lY-IFIX(Y) 

IF  (MK  NE  3)  GOTO  245  /*FOR  H  LINE 
IXO»IX  /•  SAVE  PREVIOUS  POINT 
I YO  =  I  Y 

I  PC  TR  »  0  /  +  F  0  R  FOLLOUING  POINTS 

243  CALL  I N  T  UR  T ( I Y, 6  ) 

IF  (MK  EQ  3  )  GOTO  248 
IF  (MFC  CT  17  )  N FC  =  l  7 
CALL  INTURTC NFC , 2  ) 


ML  =  3,  4  , 3 


C  PROGRAM  MAP  out 


IF  ( NFC  LT  1 >  GOTO  248 
DO  246  J  *  1 , NFC 
I=IFIX(FCODE( J)  +  0  5) 

CALL  INTUR TCI/6) 

246  CONTINUE 

248  CALL  I1TURIT(MTUNIT  ,NBLOCK  ) 

IF  (ERROR)  GOTO  989 
GOTO  200  /'♦FOR  NEXT  ITEM 
C 

260  CALL  INTURTCMR, 2  )  /'END  of  file 
CALL  NTURI TCMTUNIT /H BLOCK  ) 

IF  (ERROR)  GOTO  999 
CALL  NCLOSI 

Q *YSHO*A( ' ANY  MORE  F  I  LE  S  '  ,  1 4  ,  A  *  ND E F  ) 

IF  <  0  )  GOTO  100 
I0P  =  1 

MK  =  30 

CALL  I  NT  UR  T ( MK , 2  ) 

CALL  MTUPI T( HTUN IT , NBLOCK  ) 

IF  (ERROR)  GOTO  999 

CALL  IITUFc  (ITUNIT, ERROR,  ERROR  >  /'WRITE  FILE  MARK  TO  TAPE 
IF  (ERROR)  GOTO  999 

CALL  MTUF< MTUMI T /ERROR, ERROR  )  /*TUO  HEEDED 
IF  (ERROR)  GOTO  999 
CALL  REU(MTUNIT) 

STOP  :  123456 
C 

999  URI TE(  I ,  1090  ) 

1090  FORMAT!.  '  FATAL  MAGTAPE  ERROR') 

CALL  NCLOSI  /'CLOSE  ANY  OPEN  INPUT  FILE 

STOP  :  77777  / *  D  ON  T  REWIND  MT  -  IT  CAN  BE  INSPCTEP  FOR  DAMAGE 

END 

£»*♦••******«****»*«*•**'*****»*»»***•*♦»  *»**<  «..**♦■>»*(,*♦*•«*•»*♦'*»**'*»****» 

SUBROUTINE  MT UR  I T<  MT UN I T . NBL OC K  ) 

C  TAKES  120  CHARS  IN  IBUF,  FORMS  60  INTEGERS  AND  WRITES  THEM  TO  MAGTAPE 
LOGICAL  EOF .EOT . ERROR 
COMMOH/ERR/ERROR 
COMMON. 'INTRO C/I  BP,  IBUF<  I  ) 

C 

C 

DO  10  J  =  1 , 60 
J  B»  J ♦ J 
J A  =  JB-  1 

I  N  T  *L  S(  IBUFC  JA), 8)+RT( I 8UF ( J  6  > ,  8) 

I BUF( J  )« INT 

10  CONTINUE  /'ALL  120  ChARS  NOW  PACKED 
C 

error*  false 

CALL  II  TUR(NTUN  IT, I  BUF-60, EOT/ERROR) 

HBL0CK*NBL0CK+1 
IF  (EOT)  GOTO  900 
IF  (ERROR)  GOTO  920 
RETURN 

900  WRI TE (  1 ,  1 000  )NBL  OC  K 

1000  FORMAT!  ' MAGTAPE  EOT  AT  BLOCK’,  16) 

ERROR*  TRUE 

RETURN 

920  URI TE (  1,1001 >  H  B  L  0  C  K 

1001  FORMAT! • MAGTAPE  ERROR  IN  BLOCK’, 16) 

RETUP N 
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C  PROGRAM  MAP. OUT 


END 

C- •  —  —  •  —  *•*»»*» •*«*».,* 

C  VERSION  OF  INTURT  USED  UNTIL  NOU,  BUT  COMMENTED  OUT  NT  11-JaN-BQ 
C  SUBROUTINE  1NTWRT!I,N> 

CC  CONVERTS  AN  INTECER  TO  N  (I  TO  b)  CHARACTERS  IN  IBUF 
CC  NEGATIVE  NUMBERS  NOT  YET  AVAILABLE 
C  LOGICAL  ERROR 

C  DIMENSION  ISTRNG!  6  ) 

C  COMMON/ERR/ERROR 

C  COMMON/INTRDC/IBP,IBUF!l  > 

CC 

c  c 

C  ERROR-  FALSE. 

C  IF  (N.LT.l  OR  H  GT.6)  GOTO  990 

C  IF  <  I  LTD  GOTO  99B  /-NEGATIVE  INTEGERS  NOT  YET  AVAILABLE 

C  INT-I 

C  DO  100  J-1,N  /-PLACE  N  CHARS,  BACKWARDS,  INTO  ISTRNG 

C  1A-INT/10 

C  IR-INT-IA-10  /-REMAINDER 

C  INT-IA  /-READY  FOR  NEXT  ITERATION 

C  1STRNC! J >-IR-i 260  /-PLACE  BITS  5-8  (1011XXXX) 

C  100  CONTINUE 

C  IF  (INT.NE  0)  COTO  990  /-INTEGER  TOO  LARGE  FOR  N  CHARS 

CC 

C  DO  150  J-l.N  /-PLACE  ISTRING,  BACKWARDS,  INTO  IBUF 

C  IBUF!  IBPI-ISTRNG! N-J  +  l  > 

C  150  IBP- I  BP-1 
C  RETURN 

CC 

C  990  UR  I T  E!  1 , 1000  >1  ,  N 

C  1000  FORMAT!  ' ERROR  IN  INTURT!  '  ,216,  ’  )'  ) 

C  ERROR-TRUE 

C  RETURN 

C  END 

C-*---*»*--«-***«-*»****-*«***-***-***»---'*»--*»-***«****»*--»»**» 

C  NEU  VERSION,  TESTED  ON  VDU,  BUT  NOT  YET  ONTO  MAGTAPE 
SUBROUTINE  INTURT!  I, N  ) 

C  CONVERTS  AN  INTEGER  TO  N  !1  TO  6)  CHARACTERS  IN  IBUF 
C  SUITABLE  FOR  -3276?  <  N  <  -3276/ 

LOGICAL  ERROR 
COMMON/ERR/ERROR 
COMMON/INTRDC/IBP,  I8UF!  1  > 

C 

e 

ERROR-  FALSE 

IF  (N.LT.l  OR  N.GT  6)  GOTO  990 
NA-N  /-NO  OF  CHARS  FOR  UNSIGNED  I 
I  NT- 1 

IF  USE  0)  COTO  50  /-NOT  A  NEGATIVE  HUMBER 

NA-NA-  1 

INT--INT 

IBUF!  IBP  )-: 255  /-MINUS  SIGN 
50  IF  (HA  LT. I  )  GOTO  990 

IBP-IBP-M  /-POINT  TO  NEXT  INTWRT 

DO  100  J-l.NA  /-PLACE  NA  CHARS,  BaCKUhRDS,  INTO  IBUF 
IA-INT/1B 

IR-INT-IA-10  /-REMAINDER 

INT-IA  /-READY  FOR  NEXT  ITERATION 

IBUF!  1BP-J  )  ■  I  R  -  i  26  0  /-PLACE  BITS  5-8  tlUllXXXX) 

100  CONTINUE 


♦  *  *  *  ♦.■#*# 
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C  PROGRAM  MAP. OUT 


IP  (INT.ED  B)  RETURN 
C 

99  B  UR  I TE (  1,  1000  )1. N 
100  0  FORMAT* ' ERROR  IN  I  NT URT  (  *  .  2 1  0 .  '  )'  > 
ERROR- . TRUE . 

RETURN 

END 


i 


o  o  o 


56 


C  PROGRAM  MAP  TEX 


C  PROGRAM  MAP  TEX 

C  TO  DISPLAY  A  MAP  FILE  OH  THE  TEKTRONIX  4014 
LOGICAL  Q 
INTEGERS  DPOS 

COM MON, 'COORDS /MX , X , Y, NPTS , NPTZ, NFC, F CO  DEC  1  7  ) 
COMMOH/HEADER/I  HDUM<  100  >,  CNR*.  2  /  4  > ,  I  HDUIIA*  908  > 
fINSERT  SYSCOM>A*KEYS 


URITEC 1, 1000) 

1000  FORMAT! ' MAP  TEK  l-FEB-80') 

CALL  TKINIT 

SB  O-YSNOJAC ' FLASH  THE  SCREEN  16 , A4NDEF  ) 

IF  (0)  CALL  CLRTEK  /-FLASH  THE  SCREEN 
URITEC  1 .  1001  ) 

1001  FORMAT! ' CROSS  SIZE  (16)!') 

READ!  1 , *  )NCROSS 

C  1100  FORMAT!  1 6  ) 

CALL  NOPENI  /*  0  PEN  INPUT  FILE 
CALL  HHDIN  /*INPUT  THE  HEADER  BLOCK 
Q-Y SHOTA !' SPECIFY  CORNERS'  , 1 5,  A*DNO  > 

IF  <  NOT  Q  )  GOTO  60 

52  0-RNUM*A( ' S  U  .  X ' , 6 , A *D EC , DP  OS  ) 

IF  (  . NOT  Q  )  GOTO  52 

IF  <OPOS.LT. -1008  .OR.  DPOS  GT  5000)  GOTO  52 
SUX-FLOATC OPOS) 

53  G«RNUM*A<'$.W.  Y',6,A $t>£C,DPOS) 

IF  !  . NOT  ,Q  )  GOTO  53 

IF  (DPOS  LT  -1000  OR.  DPOS  GT  5000)  GOTO  53 
SUY-FLOAT! DPOS) 

54  Q  »R HUM 4 A! 1 N  .  E  X ' . 6 , AIDEC , DPOS  ) 

IF  <  NOT  0  )  GOTO  54 

IF  (DPOS  LE  IFIX(SUX)  .OR.  DPOS  GT  5000)  GOTO  54 
RNEX  =  FLOAT( DPOS  ) 

55  Q»RNUM*A('N  E  Y ' . 6 , A*DEC , DPOS  ) 

IF  (  NOT  0)  GOTO  55 

IF  (DPOS  LE  IFIX(SUY)  OR  DPOS  CT  5000)  COTO  55 
RHEY-FLOATC DPOS  ) 

PU-RNEX-SUX 
PH-RNE Y-SU Y 
GOTO  70 
60  SUX-0  0 
SUY-0  0 
PU-4095  0 
PH-30880 

70  U»CNR(  1. 3)-CNR<  1  ,1  )  /-UIDTH  OF  FRAME 
H-CNR( 2. 3>-CNR< 2, 1  )  /-HEIGHT  OF  FRAME 
IF  (H  LE  B  OR  U  LE  0)  COTO  950 

HSC-PU/U  / *X  SCALING  FACTOR  TO  FIT  TEKTRONIX  4014 
YSC-PH/H  / *Y  SCALING  FACTOR 
IF  (0  )  GOTO  80 

XSC-AHIN0( XSC.YSC)  /-CONVERSION  FROM  FILE  (HN>  TO  TEKTRONIX  UNITS 
YSC-KSC 

80  0-YSH04A! ' FRAME ' ,5,A*NDEF  ) 

IF  <  NOT  0  )  COTO  1 00 
00  90  J-I.5 
J  A*  J 

IF  (  J  A  EQ  5  >  J  A  =  1 
IXS-IXE  /-RUBBISH  FOR  J-I 
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C  PROGRAM  MAP . TEX 


IYS-I YE 

IXE»IFIX(CNR(  1/ JA)*XSC+0  5  +  SUK  ) 

I YE-I F 1 X C CNR< 2, J A  )*YSC*B . 5*SUY  ) 

IF  (  J  . EQ  1  )  COTO  90 
CALL  AMLCOT (29) 

CALL  POSPT C  IXS, I YS  ) 

CALL  POSPT <  IXE, I  YE  > 

90  CONTINUE 

ITEM-0  -'•COUNTS  ITEMS  (LINE/  POINT)  IH  A  SHEET 
100  CALL  NFETCH 

IF  (MX  CE.S  OR  MK.LT  0>  COTO  EBB 

IF  (MX  EQ  3  OR.  MX.EQ.4  OR  MX.EO  5)  I TEM-I TEHM 
1XP»IF1X(X*XSC*B  5+SUX) 

IYP»IFIX( Y*YSC*0  5+SUY) 

COTO  (  200.  200/ 300, 400/400  >/MK  /'•DUMMY  FOR  HK-2 
C  MX-0  OR  MX  - 1  INTERIOR  POINT  OF  A  LINE  OR  END  OF  LINE 
2BB  CALL  POSPT(  IXP,  I  YP  )  /*DR AU  TO  X,Y 
COTO  100 

C  MX  - 3  START  OF  LINE 

300  CALL  AMLCOT  (  29  )  ,‘«SET  VECTOR 

CALL  POSPT(  IXP,  I  YP  )  /»MOVE  TO  X.Y 
COTO  100 

C  MX -4  OR  S  POINT  OR  FIDUCIAL  IGNORE  FOR  THE  PRESENT 

400  CALL  PLOTCR(  IXP,  IYP,  NCR  OSS  )  ,'*PLOT  POINT  AS  A  CROSS  (MY  SUBR  ) 

COTO  100 

C  MX-E  END  OF  THE  MAP 

EBB  CALL  AML  CO  T( - 1  )  /•PURGE  BUFFER.  IE  DRAU  EVERYTHING 
URITE(  1,  1003  >MK ,  ITEM 
1003  FORMAT( ' MX-' .  IE-  '  AFTER  ITEM'.  IS) 

900  CALL  NCLOSI  /*CLOSE  INPUT  FILE 

Q-YSNO«A( ’ ANY  MORE  F I LE S ' , 1 4  -  A »NDEF > 

IF  ( Q  )  COTO  30 
STOP  : 12343E 
C 

930  URI TE(  1, 1021  ) 

1021  F ORMAT ( '  CNRS  INCORRECT.  FILE  CLOSED') 

GOTO  900 
END 

C* ••••••••»•«•••••••••••*••••••••«•••••••*••••••••••••••••»•• •••••••»••• 

SUBROUTINE  PL OT C R( IX , 1 Y , N  )  /*PLOTS  A  POINT  AS  A  CROSS.  FOR  FIPLOT 
C 
C 

CALL  AMLCOT ( 29  ) 

CALL  POSPT(  I X *N , IY  ) 

CALL  POSPT (  IX-N,  IY  ) 

CALL  AMLCOT ( 29 ) 

CALL  POSPT(  IX, I Y*N  ) 

CALL  P0SPT(  IX, I Y-N  ) 

RETURN 
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C  PROGRAM  HAP  TRhNS 


C  PROGRAM  HAP. TRANS 

C  CONVERTS  A  HAP  FILE  FROM  "IMAGE"  TO  B  N  G.  USING  A  MATRIX 
LOGICAL  Q 
INTECERM  ONUH 
DIMENSION  0CNR<2,4> 

COMHON/COORDS/HK.X.Y, HP TS, NPTJ, NFC. F CODEC  I?  > 

COMMON  ,'FILENC/NPROGC 4  I.NVERS 

COHHOH/HEADER/IHDRC  4B  >.  IPR  JC  1 ,  I  HDUH  .  SC  ALE  .  I  HDUMA(  36  ), 
•  CNR< 2.4  ). CCNR< 2,4  ),  IHDUMB( 692 ) 

COMMON/INCOH/ INDEV , INPTR,  I NS2 , I N6UF (  1 B24  ) 

•  INSERT  S YSCOM >A$KE YS 

DATA  NPROG, NVERS/BHHP. TRANS,  1/ 

C 

C 

URITEC  1,  1 0 0 B  ) 

1BBB  FORHATC ' MAP  TRANS  14-FEB-80') 

IBB  CALL  HATOPH  /»OPEN  MATRIX  FILE 
NATDEV-INDEV+4 

READ< H AIDE V, IBB  1  X  INBUFC J >, J-t , 48  > 

MR  I  TEC  1, 1001 )C INBUFC J  >, J- 1 , 48  ) 

1BB1  FORMAT <  4BA2  ) 

RE ADC  HATDEV, IBB 2  >T 1 , At  1 , A  12 
1 BB2  FORMAT C  3F2B  6  ) 

READC  HATDEV, 1BB2 >T2. A21.A22 
URI TEC  1, 1BB3)T1 . T2.A11, A12.A21 , A22 
1 BB3  FORMA TC 6F l 2 . 4  > 

CALL  NCLOSI  /*CLOSE  MATRIX  FILE 
C  CALCULATE  THE  INVERSE  OF  THE  MATRIX 
A-A 1 1  * A22  -  A2 1  * A1 2 
AIU-A22/A 
All  2- -At  2/ A 
A 12 1*- A2 1/ A 
A I22> A 1 1 /A 
C 

IIB  CALL  NOP, -HI  /•OPEN  INPUT  FILE 

CALL  HHDIN  /'•READ  THE  HEADER  BLOCK 
IF  CIPRJCl  Efl  99)  GOTO  12B  /• I HAGE  UNITS 

URITEC 1, 1BB4)IPRJC1 

1114  FORHATC 'UNSUITABLE  INPUT  HAP  PROJECTION: ' .16) 

GOTO  928 

128  URITEC 1. 1086) 

1886  FORHATC ' SUBIMACE  TOP  LEFT  CORNER') 

122  Q-RNUMfAC ' COLUMN  NUMBER 1 3 . A* DEC . DNUM  ) 

IF  ( . NOT  0  )  GOTO  122 
T1-T1-FLQATCDNUH )♦  1 . 8 
124  Q-RNUHtAC ' ROM  NUMBER ' . 1 8. A»DEC . DMUH > 

IF  C . NOT  Q  )  GOTO  124 
T2-T2-FLOATCONUH  )♦!  8 
CALL  NOPENO  /*OPEN  OUTPUT  FILE 
CALL  UPDATE  /*UPD ATE  THE  HEADER 
C  CONVERT  THE  CORNERS  AND  FIND  HAX/NIN 
XMAX-- IBBBBflfl  8 
YHAX-XNAX 
XMIN-1BBBBBBB 
YMIN-XMIN 

ROM-CNRC2,2)-CNRC2.1 )  /•HEIGHT  OF  INPUT  FILE 

DO  148  J-l. 4 

XC-CNRC 1 . J  >-T 1 

YC» R.OM-CNRC  2,  J  >-T2 

XO-XC*AIll  ♦  YC*AI 12 
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C  PROGRAM  MAP . TRANS 


YO=XC*AI 21  ♦  YC‘ A I  22 
IF  (XO  GT . XMAX)  XMAK-XQ 
IF  (XO.LT.XNIN)  XM IN* XO 
IF  <  YO  .  G  T  YHAX )  YMAX-YO 
IF  (YO.LT.YMIN)  YM 1 N* YO 
140  CONTINUE 

GCNR<  1 , 1 >=  XMI N 
GCNR( 2, 1 >*  YMI N 
CCNR* 1 , 2  )*XMIN 
GCNR* 2,2  )“ YMAX 
GCNR(  1 ,3  )»XMAX 
GCNR* 2 , 3 )• YMAX 
GCNR<  1,4  )=XMAX 
GCNR( 2,4  )* YMI N 
C 

URIT£( 1, 10SB)XHIN,YMIN,XHAX,YMAX 

1050  FORMAT*  ' GEOGRAPHICAL  CORNERS,  KM'/'S  U  .  '  , 2F 1 4  32 ' N  E .  ' , 2F1 4  3 > 
C 

CNR*  1  ,  n-0  0 
CNR*  2, 1  )  =  0  0 
CNR* 1 , 2  >*0  0 
CNR*  2 , 2 )*YNAX-YH  IN 
CNR* 1 , 3>»XHAX-XMIN 
CNR<2,3>*YHAX-YM!N 
CNR* 1 , 4  )  =  XMAX-XM IN 
CNR*  2 , 4  >*0  0 

SCALE* 1000000  0  /‘KM  STORED  AS  MM 

I  PR  JC  1  *  1  ,'*BRI  T  NAT.  GRID  ASSUMES  MATRIX  IS  IMACE  TO  GRID 

CALL  NST ART  /‘PLACE  OUTPUT  HEADER 

HK*=5 

NFC»0 

DO  160  J»1 ,4 
X  =  CNR(  1, J  ) 

Y=CNR<2. J) 

160  CALL  NST  ORE  /‘FIDUCIALS 

I  TEN* 1  /‘COUNT  INPUT  FIDUCIALS 
200  CALL  NFETCH 

IF  <HK  GE  6  >  GOTO  900 
IF  ( ITEM  CT  .4  >  GOTO  220 
ITEM*ITEHM 

IF  <  MK  EO  5)  COTO  210  /‘GOTO  200  IF  819'S  NOT  WANTED 
URI TE<  1 ,  1 005  >MK 

1005  FORMAT*' ITEM' ,13,'  HAS  Ml*-'. 16) 

GOTO  910 
210  MK°4 
MFC-1 

F CODE t  1 >-8 1 9  0  /‘STORE  OLD  FIDUCIALS  AS  8l9f 
220  XC-X-T1 

YC-R0U-Y-T2 

X«XC* A  111  ♦  YC* A 1 1 2  -  XMI H 
Y-XC*AI21  ♦  YC*A 122  -  YMI N 
CALL  NST  ORE  /‘STORE  THE  CONVERTED  CO-ORD 
GOTO  200 
C 

900  CALL  NSTORE  /‘NSTORE  THE  NK«6 
910  CALL  NCLOSO  /‘CLOSE  OUTPUT  FILE 
920  CALL  NCLOSI  /‘CLOSE  INPUT  FILE 
C 

0-YSH0$A( ' MORE  FILES' ,10, A4NDEF  ) 

IF  *0  )  GOTO  100 


U  U  u  u  u  u  *  o  o 
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C  PROGRAM  MAP  TRAMS 


STOP  < 123436 
END 

NOTE  THAT  X. V  CONVERSION  IS 
XC-X-T1 
YC-V-T2 

X«XC*AIll  ♦  YC  *  A I  12 
Y«XC*AI 21  ♦  YC • A I  22 


SUBROUTINE  MA TOP N  /*OPEN  A  MATRIX  FILE 
LOCICAL  0 

DIMENSION  I F I LE M<  2fl  > 

COMMOM/I NCOM/ INOEV 
INSERT  SYSCOM>A*KEYS 


INDEV-3 

NLEN-4B  /*UP  TO  4B  CHARS  IN  FILENAME 
1B0  0»0PNP*A< ' MATRIX  FILE'  /  11<A<READ< I F 1 LEH< NLEH. INDEV  > 
IF  (HOT .0)  GOTO  IBB 
RETURN 
END 


SUBROUTINES 


SUBROUTINE  NOPEHO  /'•TO  OPEN  A  HAP  FILE  FOR  BINARY  OUTPUT 
LOGICAL  B 

DIMENSION  I F !LEO( 2B  ) 

COMMOH/'OUTCOM/ 1 OTDE  V  <  I OT  PTR ,  I  OT  S  Z .  IOTBUF<  1824) 

INSERT  SYSCOH>A«KEYS 


I0T0EV-4  /'•OUTPUT  FILE  ON  LINE  6 

I0TS2-1B24 

NLEN-4B 

IBB  Q-OPVPSA( 'OUTPUT  FILE' <ll<AtWR!T+A*DA HF<IFILEO<NLEN< IOTDEV, 
•  AtOVAP . B  <  B )  /*PDR  3  1  Bt  P23-18 
IF  (  .  NOT . 0 >  GOTO  IBB 
RETURN 
END 
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SUBROUTINE  PERNHD  /*TO  CLEAR  HEADER  HUD  OBTHIN  PERMANENT  HEADER  DETAILS 
COMMON/HEADER/IPROJ(6 >, IDESCP  16  ), IHDUM<  1002  >  ,'*PERM  HDR 
C0MM0N/0UTC0M/I0TDUH(2).1QTS2 


DO  20  J  =  1  ,  I OTS2  /♦CLEAR  /HEADER/ 
20  I  PRO J! J  )*0 
U R  1  T E <  1  -  1000  ) 

1000  FORMAT! 'PROJECT  (12):') 

READ  (  1/1001  )(  IPROJ<J  >  .  J  »  1 . 6  ) 

1 001  FORMAT!  16A2  ) 

WRITE!  1  .  1002  ) 

1002  FORMAT! 'DESCRIPTION  (32):') 

READ  !  1/1001  )(  IDESCR! J  > .  J  =  1  .  1  6  > 
RETURN 

END 


SUBROUTINE  UPDATE  /♦TO  PLACE  DETAILS  OF  CURRENT  RUN  INTO  HEADER 
DIMENSION  IARRAY!  16  ) 

COMMON /FILEHC/NPROGV!  5  )  /♦PROG  NAME  !4>  AND  VERSION 
COMHON/HEADER/IHDR! 1024) 


DO  100  J> 1 , 360 

100  IHDR! 801-J )»1HDR! 761-J  )  /♦MOVE  DOWN  PREVIOUS  UPDATES 
DO  1 10  J»l, 3 

110  IHDR!200*J >»NPROGV!  J)  /'♦TRANSFER  PROGNAHE  AND  VERSION 
IHDR! 206>-IHDR!246>*l  /♦OUTPUT  VERS  NO 
CALL  TIMOAT!  I  ARRAY . 16)  /♦PDR  3110  P  3-54 
IHDR!207)>IARRAY( 2)  /♦DAY  OF  MONTH 
IHDR(208)»I ARRAY!  1  )  /♦MONTH  NUMBER 
IHDR! 209J-IARRAY! 3)  /♦YEAR  NUMBER  19XX 
IHDR! 210 >»I ARRAY! 4  )  /♦MINS  SINCE  MIDNIGHT 
DO  120  J«1.3 

120  IHDR!210»J )-IARRAY< 12»J>  /♦USER 
0RITE!  1 .1000  ) 

1000  FORMAT! 'DETAILS  OF  THIS  RUN  < 32 > : '  > 

READ!  1 . 1001 X IHDR! J  >. J-214, 229) 

1001  FORMAT!  1 6A2  ) 

RETURN 

END 
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SUBROUTINE  NSTART  /’•STORE  HERDER.  CLEAR  IOTBUF,  SET  IOTPTR 
COMMON /'HEADER/' IHDUH!  1024  > 

COMMON/'OUTCOM/'IOTDEV.  IOTPTR.  IOTSZ.  IOTBUF!  1  ) 

INSERT  SVSCOM>KEYS. F 


DO  IBB  J»l, IOTSZ 

IBB  IOTBUF! J >»IHDUM!  J  >  /'•TRANSFER  HEADER  TO  OUTPUT  BUFFER 

CALL  PRUF**<  KtUR I T . 10 T DE V . LOC< I OTBUF I, IOTSZ. INT L! B  ) . IRHU. ICODE ) 
IF  (  IRNU Efi . IOTSZ  AND.  ICODE. EQ.B)  GOTO  140 
WRITE! 1 ,1B8B)IRNU, ICODE 
1 BBB  FORMAT!  'PR«F»*  FAILURE  IN  NST ART  t '  ,216) 

STOP  i 111  1 
14B  IOTPTR-1 

DO  1 3B  J»l, IOTSZ 
1  SB  I OTBUF !  J  )*fl 
RETURN 
END 
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SUBROUTINE  HSTORE  /‘TO  STORE  AN  ITEM  TO  A  MAP  FILE 
DIMENSION  I R( 2  ), OUTBUF! 1  ) 

COMMON/COORDS/MK  ,  X,  Y<  NPT  S , NPT2, MFC , F CODE( 1 7 > 
COMMON/OUT COM/ IOTDEV< IOTPFS; lOTSZ, IOTBUF( 1  ) 

INSERT  SYSCON>KEYS . F 

EQUIVALENCE  (R,IR(1  >>,(IOTBUF(l  >,OUTBUF(l >> 

DATA  RMARK  / 1048576  B/ 


I SMTCH- 1  /*USED  FOR  ROUTING  THROUGH  PRUFtt 
IF  (MK  CT.6  OR.  MK  LT.B)  GOTO  900 
COTO  <  50,  900,30, 30,  50, 60  ), MK 
I T  EMSZ- 4  /* I TEMSZ  IN  INTEGERS  HK-B 

GOTO  100 

30  ITEMS2=8*NFC+NFC  /*MK»3  OR  4 
GOTO  100 

30  I T  EMSZ- 6  /*MK»1  OR  3 

GOTO  100 

60  I T EMSZ-2  /*MK-6 

100  IF  <  IOTPTR+ITEMSZ  LE .  IOTSZ+ 1  >  GOTO  2B0  /#SPACE  FOR  THIS  ITEM 
IF  (  IOTPTR  CT . lOTSZ  )  GOTO  120  /-NO  SPACE  LEFT  IN  BUFFER 
R  =  RMARK ♦?  /*“FILL”  INDICATOR 
IPTR-IOTPTR/2+l 
OUTBUF <  IPTR  >»R 

120  CALL  PRUF«*<  KSURIT. IOTOEV, LOC< IOTBUF  >, IOTS2.  IHTL<  B  >,  JRNU,  I  CODE  ) 
IF  < IRNU. EO  IOTSZ  AND.  ICODE  EQ.B)  GOTO  140 
UR  I  TE<  1 ,1000  >IRNW, ICODE 
1 B  00  FORHAT<  'PRUFtS  FAILURE  IN  HSTORE s '  , 2  I  6 > 

STOP  lllll 

140  IF  (ISUTCH  EQ.2)  RETURN 
IOTPTR-1 

DO  150  J-l. I0TS2 

130  I0tBUF(J>-8 

200  IPTR-I0TPTR/2M 

IF  (MK.EQ  0)  GOTO  300 
R-RHARK+FLOAT!  MK  > 

OU TBUF <  IPTR  )-R 
IPTR-IPTR+1 

GOTO< 5  00, 900,300, 300, 500 ,  6BB  ) , MK 

300  IRH  )-NFC+32*HPTZ 
I R( 2  )«  NPTS 
OU  TBUF (  IPTR  )«R 
IPTR-lPTR+1 

TF  (NFC  LE.0)  GOTO  500 
IF  (NFC  CT. 17)  NFC-17 
DO  320  J-l.NFC 
OUTBUF (  IPTR  )-FCODE( J  ) 

320  IPTR-IPTR+1 

300  OUT0UF (  IPTR  )-X 
IPTR-IPTR+1 
OUTBUF (  IPTR  )»Y 
IPTR-IPTR+l 
IOTPTR  - IPTR*  IPTR  - 1 
RETURN 

600  ISUTCH-2 
GOTO  120 

90B  URITE( 1 , 1081 >HK 
1001  FORMAT! ' MK- ' ,16, ' 

RETURN 
END 


IN  HSTORE  '  > 
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SUBROUTINE  NCLOSO  /*TO  CLOSE  AN  OUTPUT  NAP  FILE 
LOGICAL  0 

COHHON/OUTCON/'IOTDEV,  I  OT P Tfi , I OT S Z ,  IOTBUF<  1  ) 
INSERT  SYSCOH>Af KEYS 


0-CL0S*A<  IOTOEV)  /'•PDR  31B6  P23-19 
IF  <Q>  RETURN 
UR  1 T E<  1  ,1080  ) 

1 BBB  FORHATC 'CLOS*A  FAILURE  IN  NCLOSO  ') 
STOP  i  1  111 
END 


SUBROUTINE  NOPENI  /*TO  OPEN  A  NAP  FILE  FOR  BINARY  INPUT 
D  I  MENS  I  ON  IFILEK20) 

LOGICAL  Q 

COHNON/'INCOM/'INDEV,  INPTR,INSZ,INBUF<  1  B2  4  > 

INSERT  SYSCON>A*KEYS 


INDEV-5  /’•INPUT  FILE  ON  LINE  5 
INSZM024  /'•INPUT  BUFFER  SIZE 
NLEN-40  / *H A X  4B  CHARACTERS  IN  FILENANE 
IBB  0»OPNP*A(  '  INPUT  F I L E  '  ,  1 B  ,  A *RE AD  .  I F I L E I  .  NL E N  ,  I N D E V )  /•31BB  23-19 

IF  (  NOTH)  GOTO  IBB 
RETURN 
END 
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SUBROUTINE  NHDIN  /*READS  A  RECORD  INTO  /HEADER/ 
LOGICAL  ERROR 
CO(1HON/HEADER/INDR< 1B24> 

CONHON/INCOH/INDEV, I  HP TR , I N SZ , I N BU F(  1 B24  ) 

INSERT  SYSCOM>KEYS.F 


CALL  PRUFf»(  K»RE AD,  INDEV,LOC<  INBUF  ),  I  NS 2,  INTL<  0  ),  IRNU,  ICODE) 
IF  (ICODE  EQ  0  AND.  I RNU . E Q .  I  NS Z >  GOTO  50 
UR  I  TEC  1 , 1 020 > I RNU, ICODE 
1020  F  ORHAT<  ' PRUF44  FAILURE  IN  NHDIN)'. 2IS) 

STOP  Hill 

SB  DO  1 08  J«1,INSZ  /^TRANSFER  INTO  /HEADER/ 

100  IHDR< J )-IN8UF( J  ) 

INPTR-INSZM  /*FORCE  BUFFER  READ  AT  NFETCH 

RETURN 

END 
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SUBROUTINE  NFETCH  /*TO  FETCH  AN  ITEH  FROM  A  nAP  FILE 
DIMENSION  I  R ( 2  >/  R INBUF  < 1),ICQ0RD!42> 

COMMON/COORDS/HK.X.Y,  NPT S , NPT Z , NFC , F C ODE! 1?) 
CONMON/INCOH/INDEV, INPTRj INSZ, INBUF!  1  ) 

INSERT  SYSCOM>KE VS . F 

EQUIVALENCE  ( R , I R! 1  )  > , ( I N8UF!  1  )  , RI NBUFl  1  >  )  ,!  I  COORD C  1  ), MK  ) 
DATA  RNARK  /1048576.B/ 


DO  58  J » I >42 

58  l COORD!  J  )«B  /*CLEARS  ALL  /COORDS/ 

IF  (  INPTR  LE  INSZ  )  COTO  88 

SB  CALL  PRUF**( KtREAD, INDEV.LOC! INBUF  ),  INSZ,  INTL< B  ), IRNU, ICODE >  /•3118  3-23 
IF  < IRNU. EQ . INSZ  .AND.  ICODE. EQ.B)  GOTO  70 
WRITE! 1 , 1BBB  )IRNW /ICODE 
1 BBB  FORMAT! ' PRUF *•  FAILURE  IN  NFETCH'. 216) 

STOP  Mill 
7B  INPTRM 
8B  I PTR»I NPTR/2+1 
R*RIN8UF! IP TR ) 

IPTR-IPTR+1 

RMK-R-RMARK 

IF  (  RMK . LT . B . 9  OR.  RMK. GT. 7.1)  COTO  IBB  /*NOT  A  MK 
MK»IFIX!RNK+8. 5  > 

COTO  !  150.9  99.  130.130,  150, 16B, 6 B),MK 
IBB  MK *8 
X»R 

COTO  155 

13B  R=RINBUF( IPTR) 

IPTR»IPTR+I 

NFC-RT!  IR( I  ),5  )  /*RI GHT  5  BITS 
NPTZ=*RS!  IR<  1  ),  5)  /*H  I  CHER  PART  OF  WORD 

NPTS»IR(2  > 

IF  (NFC  LTD  COTO  150 
DO  134  J-l.NFC 
FCODE! J  )-RINBUF< IPTR) 

134  I PTR“ I P  TR ♦  1 
1  SB  X-RINBUF!  IPTR) 

IPTR=IPTR+l 
155  Y  =  RI NBUF! IPTR) 

IPTR*IPTR+1 
168  INPTR*IPTR+IPTR-1 
RETURN 

999  WRITE! 1 , 1BB1 >MK 
1BB1  FORMAT!  'MK»'  . 16.  '  IN  NFETCH') 

RETURN 

END 


SUBROUTINE  NCLOSI  /»TO  CLOSE  AN  INPUT  MAP  FILE 
LOGICAL  Q 

COMMON/INCOM/INDEV,  INPTR, I  SZ, INBUF!  1  ) 

INSERT  SYSCOM7A9KEYS 


Q»CLOS«A! INDEV  )  /*PDR  3106  P23-19 
IF  !  Q )  RETURN 
WRITE!  1 . 1 BBB  ) 

1 BBB  FORMAT ! ’ CLOStA  FAILURE  IN  NCLOSI') 
STOP  i  1111 
END 


